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1.0  Introduction 

The  ability  to  predict  accurately  the  exposure  rate  at  large 
distances  from  a  gamma  radiation  source  of  known  intensity  is  becoming 
increasingly  important  as  the  control  of  radiation  levels  in  and  around 
nuclear  facilities  becomes  more  stringent.   Of  particular  concern  are 
those  sources  of  radiation  which  result  in  a  relatively  low  exposure 
to  a  large  population  over  an  extended  length  of  time.   One  such  source 
of  low  level  exposure  is  the  reflection  of  gamma  photons  by  the  air 
above  a  radiation  source.   This  air-scattered  radiation  is  referred  to 
as  skyshine  and  may  be  an  important  consideration  in  the  design  of 
facilities  from  which  gamma  photons  may  escape  into  the  atmosphere.   For 
example,  a  reactor  containment  vessel  or  radioactive  waste  storage  building 
will  generally  have  less  shielding  in  the  roof  than  in  the  side  wails  and 
therefore,  escaping  radiation  will  have  a  preferential  upward  direction. 
In  the  design  of  such  facilities,  the  accurate  calculation  of  the  sky- 
shine  radiation  at  large  distances  from  the  facility  becomes  very  impor- 
tant. 

The  most  accurate  computational  techniques  available  today  for  the 
calculation  of  gamma  photon  transport  through  the  atmosphere  are  the 
multi-energy  group  and  multi-dimensional  transport  theory  codes.  How- 
ever, such  transport  cedes  are  very  complex  and  require  extremely 
large  computers.   Consequently,  the  use  of  these  photon  transport  codes 
for  skyshine  calculations  is  prohibitively  expensive  for  routine  analysis, 
often  requiring  several  thousand  dollars  worth  of  computer  time  for  a 


single  run  if  a  high  level  of  accuracy  is  required. 

Because  of  the  large  financial  burden  associated  with  the  routine 
use  of  these  state-of-the-art  transport  codes,  much  effort  has  been 
expended  in  the  development  of  simplified,  approximate  models  to  be 
used  specifically  in  the  calculation  of  skyshine  exposure  rates.   The 
most  popular  of  these  simplified  models  are  based  on  the  point  kernel 
technique  which  calculates  a  first  collision  source  distribution  in  the 
atmosphere  around  the  source.   The  angular  distributions  of  the  secondary 
photons  are  given  by  the  Klein-Nishina  cross-section.   Buildup  factors 
are  then  used  to  account  for  subsequent  scattering  between  the  location 
of  the  first  interaction  and  the  location  at  which  the  exposure  rate  is 
to  be  calculated. 

Although  such  simplified  calculational  models  are  quite  inexpensive 
to  use  compared  to  the  sophisticated  transport  codes,  the  accuracy  of  these 
models,  particularly  for  large  source/detector  distances,  is  generally  un- 
known.  A  working  committee  was  recently  formed  by  the  Nuclear  Safety 
Research  Association  of  Japan  (NSRA)  because  of  renewed  interest  in  the  con- 
trol of  low  level  exposure  to  the  general  population  .   The  membership  of 
this  committee  consisted  of  representatives  from  various  Japanese  utilities, 
consulting  firms  and  research  organizations.   Skyshine  exposure  rate 
predictions  for  various  point  source  configurations  were  calculated  by 
members  of  this  committee  using  several  different  point  kernel  schemes. 
The  results  obtained  by  this  committee  exhibited  large  discrepancies 
between  the  values  calculated  by  the  different  schemes.   Often  a  factor 


of  almost  seventy  would  separate  the  extreme  results  at  the  same 
source/detector  distance. 

The  large  disagreement  between  many  of  these  results  obtained  by 
the  working  committee  indicated  a  large  uncertainty  in  the  ability  to 
predict  skyshine  exposure  rates.  This  uncertainty  and  the  importance 
of  the  skyshine  problem  prompted  the  committee  to  design  and  sponsor 
a  benchmark  skyshine  experiment  to  provide  definitive  data  to  compare 
with  calculated  predictions.  Although  other  skyshine  experiments  have 
been  previously  performed,  the  source-detector  geometries  employed  in 
these  studies  were  generally  so  complex  that  three-dimensional  skyshine 

models  would  be  required  for  an  accurate  description  of  the  experimental 

1.2 
configurations 

A  benchmark  skyshine  experiment,  designed  and  sponsored  by  the 
NSRA,  was  performed  from  May  through  November  1977,  at  the  Kansas  State 
University  Nuclear  Engineering  Shielding  Facility  (KSUNESF)  which  is 
located  approximately  five  miles  west  of  the  Kansas  State  University 
campus  in  Manhattan,  Kansas.  A  map  of  the  shielding  facility,  which 
is  approximately  175  acres  in  extent,  is  shown  in  Fig.  1. 

To  obtain  skyshine  radiation  fields  of  sufficient  intensity  to  be 
measured  over  the  700  meter  measurement  baseline,  three  Co-60  sources 
with  nominal  strengths  of  10,  250  and  5000  Ci  were  employed  in  specially 
fabricated  skyshine  irradiators.  An  approximately  cylindrical  concrete 
silo  was  built  with  walls  sufficiently  thick  so  that  the  direct  pene- 
tration component  of  the  exposed  sources  would  be  negligible.  Measure- 


ments  of  skyshine  radiation  fields  from  three  different  source  configu- 
rations were  made,  first  with  the  source  radiation  collimated  into  a  150° 
upward  directed  conical  beam  and  later  with  the  sources  shielded  by  8  and 
16  inches  of  concrete. 

Spectral  measurements  of  the  skyshine  radiation  field  were  taken 
with  a  Nal(Tl)  scintillation  spectrometer.  The  Nal  crystal  was  housed 
in  an  eight  ton, lead,  steel  and  concrete  collimator  assembly  to  reduce 
the  background  count  rate  and  define  a  collimation  angle  for  the  sky- 
shine radiation.  The  collimator  assembly  along  with  the  multi-channel 
analyzer  and  supporting  electronics  were  contained  in  a  weathertight 
metal  shed  and  mounted  on  a  twelve  ton  capacity,  8'  by  29'  semitrailer. 
The  trailer  and  detection  assembly  were  moved  over  the  700  meter  measure- 
ment baseline  by  a  five  ton  capacity  tractor-truck.   Additionally,  total 
exposure  rate  measurements  were  made  with  a  high  pressure  ionization 
chamber  at  source/detector  distances  ranging  from  30  to  700  meters. 

This  work  reviews  and  examines  the  experimental  facilities,  the 
experimental  methods,  and  the  original  analysis  of  the  benchmark  data 

which  was  performed  by  Radiation  Research  Associates,  Inc.,  of  Fort 

3 
Worth,  Texas  .   Due  to  the  large  magnitude  of  the  necessary  facilities, 

much  effort  (including  the  author's)  was  devoted  to  building  the  source 

silo  and  the  detector  collimator  assembly.  The  author  also  played  a 

major  role  in  the  integration,  calibration,  and  implementation  of  the 

gamma  spectrometer.   These  activities  are  described  in  Chapter  2. 


An  independent  analysis  of  the  skyshine  spectral  measurements 

is  presented  in  Chapter  3.   The  measured  spectra  were  unfolded  with 

4 
the  SEGO  unfolding  code  and  compared  to  the  results  of  detailed  trans- 
port theory  calculations  obtained  with  DOT-3.5  two-dimensional 
transport  code.  An  estimate  of  the  accuracy  of  the  transport  code's 
predictive  abilities  are  obtained  through  this  comparison. 

The  high  pressure  ionization  chamber  measurements  must  be  corrected 
for  the  non-ideal  energy  sensitivity  of  the  detector.   Methods  are 
presented  for  obtaining  these  corrections  and  the  corrected  measurements 
are  then  compared  to  predictions  calculated  from  transport  theory. 
These  results  and  comparisons  are  presented  in  Chapter  4.   Additionally, 
the  unfolded  Nal  spectra  are  used  for  the  first  time  to  obtain  an  inde- 
pendent determination  of  the  total  skyshine  exposure  rate.  These  values 
are  also  compared  to  the  high  pressure  ionization  chamber  measurements 
in  Chapter  4. 


2.0  Experimental  Facilities  and  Methods 

A  great  deal  of  effort  was  expended  in  the  construction  of  the 
facilities  necessary  for  this  experiment.   Because  of  the  large  distances 
over  which  skyshine  measurement  were  to  be  taken,  a  wide  range  of  source 
strengths  were  used.  This  necessitated  a  large  irradiation  unit  to  pro- 
vide sufficient  shielding  of  the  sources  and  to  accommodate  the  different 
source  configurations  utilized.   The  preparatory  activities  centered  around 
three  main  projects: 

1)  the  integration  and  calibration  of  the  spectrometer, 

2)  the  construction  of  the  detector  cavity  and  collimator,  and 

3)  the  construction  of  the  source  silo  and  its  associated  components. 
This  chapter  will  describe  the  above  items  and  their  fabrication. 

2.1  Gamma-Ray  Spectroscopy  System 

Due  to  the  relatively  large  source  detector  distances  over  which  sky- 
shine  measurements  were  to  be  made,  the  foreground  to  background  signal 
ratio  was  expected  to  approach  unity  at  the  remote  measurement  locations. 
Because  of  this  consideration,  a  Nal(Tl)  gamma-detection  system  was 
chosen  over  a  semiconductor.   Although  the  resolution  of  a  semiconductor 
based  spectrometer  generally  exceeds  that  of  a  Nal  spectrometer, 
the  comparatively  low  efficiency  of  the  former  system  precluded  its 
use  in  measuring  the  anticipated  low  intensity  radiation  spectra. 
Supporting  the  five  inch  diameter  by  four  inch  deep  cylindrical  Nal  crystal 
was  an  RCA  model  4525  photomultiplier  tube  (PMT)  plus  an  ORTEC  model  276 


preamplifier,  amplifier  (model  450)  and  high  voltage  supply  (model  456)  . 
Additionally,  a  Canberra  model  8180  multichannel  analyzer  (MCA)  and  a 
model  8531A  tape  control  unit  with  a  PERTEC  model  7820-9  tape  drive 
were  used  to  collect  and  store  the  data.  An  ELGAR  model  3000  line  con- 
ditioner was  used  to  supply  a  constant  line  voltage  to  the  necessary 
components.  The  entire  system  is  shown  schematically  in  Fig.  2. 

Although  the  Canberra  MCA  had  a  built-in  amplifier  and  pile-up 
rejector  circuit,  these  components  were  designed  primarily  for  a  GeLi 
based  spectrometer  system  and  did  not  perform  adequately  with  Nal 
pulses  because  of  the  very  different  pulse  shaping  requirements  for  the 
two  detectors.   Consequently,  an  external  amplifier  with  a  variable 
pulse  shaping  feature  was  used  and  the  internal  amplifier  of  the  MCA  was 
bypassed.  With  the  external  amplifier  a  wide  variety  of  time  constants 
were  available  as  well  as  a  choice  between  unipolar  and  bipolar  output 
pulses.  All  of  the  above  parameters  as  well  as  the  high  voltage  supplied 
to  the  PMT  needed  to  be  selected  to  accommodate  a  wide  range  of  photon 
energies  and  count  rates. 

2.1.1  Choice  of  Spectrometer  System  Parameters 
The  system  operating  parameters  were  chosen  so  that  the  Nal  reso- 
lution would  not  be  degraded  while  the  system  was  experiencing  a  high 
count  rate  situation.  With  bipolar  pulses  under  high  count  rate  con- 
ditions both  sides  of  a  photopeak  may  become  broadened  as  a  result  of 
pulse  superposition,  thus  causing  the  resolution  of  the  system  to  de- 
crease. Alternatively,  with  unipolar  pulses  only  the  high  energy  side 


9 

tu 
> 

•H 
U 
Q 

t 

o 

CJ  CS 

0  as 

■U   P-. 

h 

a    • 

a.  tj 

o 

s 

u 
c 

c 

1- 
4J 

c 
c 
c 

> 

o 

0) 

a 

(0 

CO 
H 

•r- 

e 

of 

O     N 

M    = 

«j 

H 

0) 
3 

--t 

0 

< 

!     O 

3 

u      at 

<u 

— 

o 

H 

-C 

N      o 

t 

i 

^    S 

■H 

•H 

c      u 
<     u 

1      <n 

0 

u 

< 

• 

2              X> 

4J 

G 

1 

0 

u 

0 

o 

^3   tH 

C     11 

4) 

a 

u 

d) 

ct) 

a. 

en 

>> 

w 

o 

u 

V 

in 

•H 

u  <■ 

N 

tH 

H 

H   — 1 

OS   u 

O   T3 

m 

i 

a 

*-Q 

M 

1 

£ 

0 

> 

in 

1 

—t 

M 

00 

•H 

a 

ON 

M 

U 

0 

1 

H   H 

as   i) 

0    "3 

S3  ^ 
H 

9) 

35 

cu 

M 

O   "3 

3 

CU 

0 

^ 

oi 

fa 

-a 

^c 

3 

m 

4J 

in 

H 

<r 

rH 

U   <T 

3 

0 

> 

fa3 

si    u 

as 

« 

w    0 

<T       M 

X       fl 

:        -z. 

m 

10 


of  the  peak  is  broadened  and  is  therefore  more  easily  discernible  by 
visual  inspection.   Furthermore,  with  unipolar  pulses  one  has  the  choice 
of  using  a  "base  line  restore"  option  in  the  amplifier.  This  option 
forces  the  trailing  edge  of  the  pulse  to  zero  thus  reducing  the  proba- 
bility of  a  second  pulse  being  superimposed  on  top  of  the  trailing  edge 
of  the  previous  pulse  and  thereby  giving  too  high  an  amplitude  for  the 
second  pulse.  For  the  high  count  rates  expected  in  this  experiment,  the 
use  of  a  pile  up  rejection  circuit  was  most  desirable. 

The  magnitude  of  the  voltage  applied  to  the  PMT  also  affects  the 
resolution  of  the  system.   In  general,  higher  resolution  is  achieved 
with  higher  voltages  although  under  high  count  rate  conditions  too  high 
a  PMT  voltage  will  cause  a  voltage  drop  in  the  last  dynode  stages  (be- 
cause of  the  high  current  flow)  which  in  turn  causes  a  downward  shift 
in  the  photopeak  position  towards  lower  energy.  The  900  V  PMT  voltage 
chosen  for  this  system  was  found  not  to  degrade  the  spectrometer  resolution 

significantly,  while  at  the  same  time,  spectral  shifting  was  less  than 

4 
1  channel  for  count  rates  up  to  3.0x10  counts  per  second  (cps). 

The  time  constant  of  the  amplifier,  which  regulates  the  width  and 

rise  time  of  the  amplifier  output  pulse,  is  an  important  parameter 

which  can  greatly  affect  the  resolution  of  the  spectrometer.   It  was 

determined  experimentally  that  in  general,  the  resolution  of  the  system 

deteriorates  as  the  value  of  the  time  constant  increases.  Also  associated 

with  the  larger  time  constants  is  a  significant  increase  in  the  dead  time 

of  the  MCA.   However,  if  too  short  a  time  constant  is  chosen  the  resolu- 
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tion  of  the  system  is  adversely  affected.   It  was  observed  that  slightly 
more  detail  was  discernible  from  spectra  obtained  with  a  0.25  usee 
amplifier  time  constant  than  with  other  values,  and  it  was  this  value  that 
was  chosen  for  normal  operating  conditions. 

2.1.2  Spectrometer  Calibration 

For  the  resulting  energy  spectrum  obtained  by  the  MCA  to  be  useful, 
the  energy  corresponding  to  any  channel  must  be  accurately  known.   For 
the  present  application  an  energy  bin,  or  channel,  width  of  5  keV  was 
specified  for  all  spectra,  which  were  to  be  512  channels  in  length. 
Although  a  variety  of  combinations  of  the  amplifier  gain  and  MCA  con- 
version gain  will  yield  this  desired  result,  the  input  pulses  to  the 
analyzer  should  be  as  small  as  possible  to  minimize  dead  time  while 
using  an  MCA  conversion  gain  value  sufficiently  large  to  obtain  the 
desired  spectral  spread. 

Since  the  gain  of  a  scintillation  system  is  affected  by  temperature, 
minor  adjustments  were  necessary  in  the  field.   The  amplifier  fine  gain 
adjustment  was  used  to  correct  for  these  small  temperature  induced  gain 
shifts.   The  gain  shifts  experienced  in  the  field  were  generally  less 
than  one  channel  at  660  keV. 

The  preliminary  calibration  was  performed  using  Am-241  (59.54  keV) , 
Ba-133  (SO  keV)  and  Cs-137  (661.6  keV)  sources.   These  sources  were  used 
to  determine  the  amplifier  gain  and  the  ADC  baseline  setting  (which 
fixes  the  channel  corresponding  to  zero  energy).   After  the  initial 


12 


calibration  a  mixed  isotope  source  was  used  to  confirm  the  linearity 
of  the  calibration  over  a  large  energy  range.   This  energy  range 
extended  from  59.54  keV  to  2.505  MeV.   The  result  of  this  calibration 
test  is  shown  in  Fig.  3  where  only  slight  deviations  from  the  ideal 
linear  model  can  be  seen  below  400  keV  while  above  this  point,  the 
response  is  extremely  linear. 

2.1.3  System  Response  as  a  Function  of  Count  Rate 
In  order  to  determine  the  maximum  counting  rate  which  could  be 
tolerated  by  this  system  without  adversely  affecting  its  calibration 
or  linearity,  the  response  of  the  spectrometer  to  incident  fluxes  of 
varying  intensities  was  studied  in  great  detail.  A  nominal  3.8  mCi 
Cs-137  source  was  positioned  at  various  distances  from  the  detector 
face  ranging  from  122  to  1036  cm.  A  precision  pulser,  (ORTEC  model 
419)  shown  in  Fig.  2,  was  also  used  to  test  the  accuracy  of  the  internal 
clock  and  dead  time  correction  circuits  of  the  analyzer.   The  pulser 
gain  was  adjusted  so  that  the  pulser  peak  fell  well  above  the  Cs-137 
photopeak.   The  area  of  the  pulser  peak  is  then  directly  related  to  the 
true  counting  time. 

Standard  512  channel  spectra  were  taken  at  regular  source-detector 
distances.  A  background  spectrum  was  also  measured  curing  the  experiment 
and  was  subtracted  in  a  channel-by-channel  fashion  from  each  of  the 
Cs-137  spectra  after  appropriate  normalization  to  correct  for  different 
counting  times.  The  background  corrected  Cs-137  photopeaks  were  then 
fitted  in  a  least-squares  manner  to  a  Gaussian  curve.   This  Gaussion  fit 
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yielded  values  for  the  photopeak  centroid  as  well  as  the  area  of 
the  fitted  curve.  Also  evaluated  for  each  spectra  were  the  area  of 
the  total  spectrum  (neglecting  the  pulser  peak)  from  which  the  total 
count  rate  is  calculated,  and  the  area  of  the  pulser  peak  (including 
any  pile-up  tail)  from  which  the  true  counting  time  may  be  found. 

In  Fig.  4  the  centroid  of  the  photopeak  is  shown  as  a  function  of 

the  total  counting  rate.  As  can  be  seen,  no  significant  gain  shift 

4 
occurred  for  total  count  rates  below  3.0x10  cps. 

Another  parameter  of  interest,  the  full-width-at-half -maximum 

(FWHM)  is  calculated  from  the  least  squares  fitted  Gaussian  to  the 

photopeak.  For  an  ideal  spectrometer,  the  resolution  (or  FWHM)  should 

be  independent  of  the  count  rate.   The  results  for  the  present  system 

are  shown  in  Fig.  5  as  a  function  of  the  total  counting  rate.  Again,  no 

significant  degradation  in  the  resolution  is  observed  below  count  rates 

4 
of  3.0x10  cps. 

Another  method  used  to  test  for  photopeak  shifting  under  high  count 

rate  situations  is  to  take  a  channel-by-channel  ratio  of  the  photopeaks, 

choosing  one  spectrum  obtained  with  a  low  count  rate  as  a  reference.  The 

results  of  this  method  of  comparison  are  shown  in  Fig.  6  where  a  nonzero 

slope  indicates  a  shift  in  the  photopeak  position.   No  significant  shifting 

4 
is  observed  below  total  count  rates  of  3.0  x  10  cps. 

Analysis  of  the  pulser  peak  and  its  pile-up  tail  indicated  that  the 

error  in  the  dead  time  correction  was  less  than  2%  in  15  out  of  16  (94%) 

4 
tests  below  3.0x10  cps.   It  was  also  verified  that  there  was  a  definite 
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linear  relationship  between  dead  time  and  the  total  count  rate. 

The  various  values  of  the  operating  parameters  used  in  this 
experiment  are  listed  below: 

1)  Unipolar  amplifier  pulses  with  0.25  usee  time  constant 

2)  Base  line  restore  (BLR)  set  to  HI 

3)  Amplifier  coarse  gain  at  100 

4)  Amplifier  fine  gain  at  11.22  (variable) 

5)  PMT  high  voltage  at  900  volts 

6)  MCA  conversion  gain  at  1024/2 

7)  ADC  baseline  set  at  4.96. 

With  these  settings,  stable  system  response  can  be  expected  from 
30  KeV  to  2.6  MeV  and  for  total  count  rates  less  than  30,000  cps. 

2.2  Detector  Collimator  Assembly 

One  of  the  major  construction  activities  associated  with  this 
experiment  was  the  design  and  fabrication  of  a  well  shielded  detector 
cavity  and  collimator  or  shadow  shield  .   Although  the  reduction  of  the 
detector  geometry  from  4ir  to  the  collimated  configuration  complicated 
the  prediction  and  analysis  of  the  results,  this  feature  was  necessary 
to  increase  the  signal  to  background  ratio  in  order  to  extend  the  base- 
line over  which  skyshine  measurements  could  be  taken.  Moreover  the 
collimated  skyshine  spectra  were  thought  to  provide  a  more  severe 
test  than  4ir  spectra  would  have  provided  for  the  transport  calculational 
models  whose  predictive  abilities  were  the  overall  objective  of  this 
project. 

To  eliminate  the  need  for  long  signal  cables,  the  detector  cavity/ 
collimator  and  MCA  systems  were  mounted  on  an  8'  x  29'  12  ton  capacity 
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flat  bed  taller  to  be  pulled  by  a  5  ton  tractor  truck.   Since  the  esti- 
mated weight  of  the  collimator  assembly  was  well  over  8  tons,  equipment 
of  this  magnitude  was  necessary.   Integration  of  the  two  systems  in 
this  manner  provided  the  further  advantage  of  requiring  only  one  temper- 
ature control  system  to  ensure  the  response  stability  of  the  spectrometer. 

2.2.1  Construction  of  Collimator  Assembly 
The  collimator  assembly  consisted  of  two  main  chambers,  the 
details  of  which  are  shown  in  Fig.  7.  This  assembly  was  constructed 
on  the  rear  of  a  flat  bed  trailer  with  the  cavity  axis  parallel  to 
the  trailer  axis  and  the  collimator  opening  oriented  towards  the  rear 
of  the  tailer.  The  inner  chamber  or  detector  cavity  was  a  24"  cube  in 
which  the  Nal  crystal  was  to  be  positioned.   The  outer  chamber,  to 
which  the  shadow  shield  was  attached,  was  3'  square  by  20"  deep.  The 
walls  of  both  chambers  were  constructed  of  5/16"  steel  plate  and  angle 
iron.  The  shadow  shield  was  composed  of  standard  sized  lead  bricks 
(2"  x  4"  x  8")  in  an  angle  iron  frame  as  shown  in  Plate  1.  An  8" 
square  access  port  penetrated  the  rear  face  of  the  detector  chamber  to 
be  used  primarily  as  a  cable  conduit  but  was  large  enough  to  permit 
removal  of  the  entire  Sal  crystal/PMT/preamplif ier  assembly  if  necessary. 
The  detector  cavity/collimator  assembly  is  shown  in  Plate  2. 

In  an  effort  to  eliminate  any  radiation  entering  the  detector 
cavity  except  through  the  shadow  shield  aperture,  nominal  8"  lead 
slabs  were  placed  on  the  top  and  bottom  of  the  inner  chamber  while 
standard  lead  bricks  were  stacked  around  the  sides  and  rear  as  shown 
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in  Plate  3.   In  addition  to  the  lead  shielding,  concrete  blocks  were 
stacked  around  the  sides  and  top  of  the  cavity  to  a  minimum  thickness 
of  eight  inches.   In  order  to  hold  the  shielding  materials  in  place 
during  the  course  of  the  measurements,  plywood  panels  were  used  to  en- 
case the  sides  of  the  assembly  which  were  in  turn  held  in  place  by 
circumferential  cables. 

The  detector  cavity/collimator  assembly  not  only  served  to  reduce 
the  background  radiation  levels  by  a  factor  of  five,  but  also  proved 
useful  in  providing  a  large  thermal  inertia  for  the  crystal/PMT/preamp 
assembly.   Thus  in  the  case  of  a  failure  of  the  environmental  control 
system,  the  crystal  assembly  would  be  well  insulated. 

2.2.2  Nal  Crystal  Positioning 

Styrofoam  was  chosen  to  support  the  detector  assembly  inside  the 
detector  chamber  because  of  its  good  thermal  insulating  properties, 
low  density,  and  ability  to  isolate  the  crystal  from  minor  vibrations. 
The  crystal  assembly  was  placed  in  a  cradle  made  out  of  3"  thick 
styrofoam  slabs,  the  bottom  half  of  which  is  shown  in  Plate  4.   Small 
styrofoam  wedges  were  used  to  firmly  fix  the  detector's  position  with 
respect  to  the  shadow  shield  aperture.   Finally,  to  protect  the  front 
face  of  the  detector,  a  two-inch  thick  styrofoam  cover  which  had  been 
hollowed  out  to  a  thickness  of  less  than  one  inch  in  the  central  portion 
was  placed  in  the  aperture  of  detector  chamber.   Thus,  the  detector 
system  was  completely  insulated  without  significant  attenuation  of  an 
incident  radiation  field. 
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2.2.3  Environmental  Control 

Since  the  response  of  the  spectrometer  was  affected  by  temperature, 
environmental  control  of  the  crystal  assembly  as  well  as  the  analyzer 
system  was  imperative.   To  accomplish  this,  two  metal  utility  sheds 
were  joined  together  and  attached  to  the  trailer  in  such  a  way  as  to 
cover  the  entire  detector  cavity/ collimator  assembly.  There  was  also 
sufficient  room  behind  the  detector  chamber  for  the  supporting  data 
acquisition  instrumentation  and  its  operating  personnel.  Large  wooden 
doors  were  installed  in  the  rear  end  of  the  shed  so  that  when  opened, 
the  aperture  of  the  shadow  shield  would  be  unobstructed. 

Cooling  for  the  instrumentation  and  detector  assembly  was  provided 
by  two  thermostatically  controlled  air  conditioning  units  with  a  total 
capacity  of  10,000  Btu/hr.   These  were  installed  in  the  walls  of  the 
environmental  shed  as  shown  in  Plate  5 .   To  reduce  the  amount  of  heat 
flow  into  the  shed,  styrofoam  panels  were  attached  to  the  walls  and  roofs, 
all  joints  were  caulked  and  weatherstripping  was  used  around  the  entrance 
door.  Finally  a  2"  thick  styrofoam  panel  was  placed  in  the  shadow  shield 
opening  and  plastic  sheeting  was  stretched  across  the  collimator  door 
opening  thereby  completely  sealing  the  rear  portion  of  the  shed.   These 
features  can  be  seen  in  Plate  5.  With  these  measures,  a  constant 
temperature  could  be  maintained  on  all  but  the  hottest  days.   The  shed 
besides  providing  temperature  control  also  provided  a  dry  and  dust-tree 
environment  for  the  spectroscopy  system. 
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Plate  5. 


Rear  of  detector  trailer  showing  collimator  opening 
and  trailer  leveling  system. 
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The  complete  tractor-trailer  detection  assembly  is  shown  in 
Plate  6. 

2.2.4  Detector  Cavity/Collimator  Performance 

To  test  the  performance  of  the  collimator  assembly,  several  energy 
spectra,  or  response  functions,  were  taken  using  four  National  Bureau  of 
Standards  (NBS)  calibration  sources  (Am-241,  Hg-203,  Cs-137,  and  Co-60) . 
Figure  8  shows  a  Co-60  response  function  with  the  source  located  on 
the  detector  axis  69.85  cm  from  the  crystal  face.   Although  the  high 
energy  portion  (>.8  MeV)  of  the  spectrum  appears  normal,  one  immediately 
notices  several  large  scattering  peaks  below  400  KeV.   These  features 
are  attributable  to  photons  which  have  been  multiplely  scattered  from 
the  steel  walls  of  the  detector  cavity.   Such  low  energy  structures 
were  highly  undesirable  since  the  skyshine  spectra  were  expected  to 
be  composed  of  predominantly  low  energy  photons. 

To  minimize  the  scattering  within  the  detector  collimator/shield 
assembly,  two  modifications  to  this  assembly  were  made.   First  the 
walls  of  the  detector  chamber  and  the  front  portion  of  the  outer 
chamber  were  lined  with  1/2"  lead  sheeting  in  order  to  reduce  the 
gamma-ray  reflection  from  the  walls.   Secondly  an  inner  lead  collimator 
with  a  12"  square  opening  was  placed  in  front  of  the  aperture  of  the 
detector  cavity  in  order  to  reduce  the  number  of  photons  admitted  to 
the  detector  cavity.   This  apparatus  is  shown  in  Plate  7.   Also  during 
these  alterations,  the  crystal  assembly  was  moved  forward  in  the  de- 
tector chamber  to  increase  the  collimation  angle.   The  final  detector 
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cavity/collimator  is  shown  in  Fig.  9.   Figure  10  shows  a  Co-60  response 
function  taken  with  the  new  collimator/shield  assembly  superimposed  on 
Fig.  8  where  the  affect  of  the  new  collimator  design  and  crystal  position 
is  apparent.   The  previous  collimator  scatter  peaks  are  no  longer  present 
and  the  Compton  plateau  more  nearly  resembles  the  ideal  Nal  response. 

To  determine  the  effectiveness  of  the  collimator  in  defining  the 
detector  geometry,  a  Cs-137  calibration  source  was  moved  vertically 
and  horizontally  away  from  the  cavity  axis  parallel  to  the  plane  of  the 
collimator.  At  various  radial  distances  from  the  detector  axis,  response 
functions  were  measured.   Ideally,  as  the  source  moves  across  the 
collimator  opening  neither  the  total  count  rate  nor  the  spectral  shape 
should  change.   However,  as  the  source  moves  into  positions  where  part 
or  all  of  the  crystal  is  blocked  by  the  shadow  shield  (see  Fig.  12), 
the  total  count  rate  should  decrease  dramatically  and  eventually  approach 
that  of  background. 

Some  small  deviations  from  the  ideal  behavior  were  noted  and  can 
be  seen  from  the  collimator  traverse  results  shown  in  Fig.  11.  Since 
the  crystal  size  was  relatively  large  compared  to  the  source-detector 
distance,  various  portions  of  the  crystal  are  shielded  by  the  shadow 
shield  at  different  positions  as  shown  in  Fig.  12.  Consequently,  the 
total  count  rate  does  not  drop  off  as  sharply  as  would  be  experienced 
with  a  point  detector.  A  second  non-ideal  response  contribution  is 
that  as  the  sources  moves  away  from  the  detector  axis  a  greater  cross 
sectional  area  of  the  crystal  is  exposed  causing  the  total  count  rate 
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Fig.  11.  Results  of  vertical  and  horizontal  traverses 
of  the  collimator  opening. 
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to  rise  slightly.  A  correction  factor  for  this  phenomenon  will  be 
discussed  in  Chapter  4. 

2.3  Source  Silo  and  Associated  Components 

The  second  major  construction  endeavor  undertaken  for  this  pro- 
ject was  the  fabrication  of  the  source  silo  and  its  supporting  com- 
ponents  .  The  silo  was  designed  to  house  the  two  Co-60  source  irradiators 
used  in  the  project  and  to  provide  sufficient  shielding  so  that  the  com- 
ponent of  radiation  penetrating  the  walls  would  be  negligible  compared 
to  the  skyshine  component.   In  addition  to  the  silo,  collimator  wedges 
were  fabricated  and  used  to  define  the  source  collimation  angle  in  the 
unshielded  skyshine  portion  of  the  experiment.   Concrete  roof  slabs 
were  fabricated  and  employed  in  the  two  shielded  measurement  con- 
figurations. 

2.3.1  Source  Silo 

To  facilitate  the  numerical  modeling  of  the  source  configuration, 
a  silo  with  approximately  an  annular  cross  section  was  desired. 
Practical  considerations  dictated  the  use  of  concrete  forms  with  flat 
sides  resulting  in  a  silo  in  the  shape  of  dodecagon  as  shown  in  Fig.  13. 
To  support  the  mass  of  the  silo  and  the  source  casks,  a  9"  thick  base 
pad  was  poured  containing  large  amounts  of  reinforcing  steel  and  wire 
mesh.  On  the  four  sides  of  the  base  pad,  4"  thick  slabs  were  poured 
to  provide  easy  access  to  the  silo.   The  completed  silo  is  shown  in 
Plate  8. 


■ 
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p"   MANHOLE  ENTRANCE 


Fig.  13.   Cross  sectional  view  of  the  source  silo 
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Plate 


Completed  source  silo.   A-frarae  crane  is  used  for 
positioning  of  roof  shields. 
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The  top  and  bottom  halves  of  the  silo  wall  were  poured  in  two 
distinct  stages  using  step  joints  between  both  halves  as  well  as  at 
the  wall  base.  The  concrete  was  also  thoroughly  vibrated  to  reduce  the 
possibility  of  voids  within  the  wall.   Several  small  penetrations  at 
the  bottom  of  the  silo  wall  were  made  to  provide  access  ports  for 
electrical  cables  and  positioning  of  the  source  irradiators. 

A  "line  of  sight"  penetration  was  made  in  the  top  front  side  of 
the  wall  so  that  the  Co-60  sources  could  be  viewed  directly  from  some 
of  the  closer  measurement  positions.   Normally  this  opening  was  com- 
pletely filled  with  standard  sized  concrete  blocks  which  were  removed 
only  for  the  source  calibration  measurements. 

2.3.2   Collimator  Wedges 

For  the  unshielded  skyshine  measurements  it  was  required  that 
the  radiation  field  be  collimated  into  an  upward  150°  conical  beam. 
This  was  accomplished  by  positioning  wedge  shaped  concrete  blocks 
around  the  top  of  the  source  silo  wall  thus  forming  the  perimeter  of 
a  cone  with  a  150°  angle  whose  projected  apex  coincided  with  the  position 
of  the  raised  Co-60  sources.   In  this  way  the  collimated  beam  strength 
could  be  changed  by  simply  moving  the  desired  source  cask  to  the  center 
position  of  the  silo  so  that  the  raised  source  would  be  at  the  apex  of 
the  cone.   With  this  technique,  no  adjustment  of  the  collimator  wedges 
would  be  required. 

The  dimensions  of  the  collimator  wedges,  shown  in  Fig.  14,  were 
chosen  so  that  48  of  them  would  form  an  annulus  with  an  8 '2"  inside 
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24" 


Fig.  14.   Dimensions  of  collimator  wedges. 
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diameter.  Arranged  in  this  manner,  a  colliniation  cone  with  a  150°  angle 
would  be  obtained.   To  minimize  radiation  streaming  between  the  wedges 
and  direct  penetration,  concrete  and  lead  bricks  were  placed  around 
the  outside  of  the  wedges  in  such  a  way  as  to  not  affect  the  collimated 
radiation  field.   The  final  collimation  angle  was  measured  and  found  to 
be  within  0.2°  of  the  design  value. 

2.3.3  Silo  Roof  Shields 

The  shielded  portion  of  the  experiment  called  for  the  radiation 
sources  to  be  covered  by  20  and  40  cm  thick  shields  of  concrete.   Con- 
sequently 4'  x  12'  x  4"  concrete  slab  shields  were  constructed  by 
pouring  concrete  into  steel  frames  with  stepped  sides  thus  allowing  adja- 
cent slabs  to  overlap.  To  provide  the  needed  mechanical  strength  so  that 
the  slabs  could  withstand  being  moved  repeatedly,  minimal  amounts  of 
structural  steel  and  wire  mesh  were  imbedded  in  the  concrete. 

To  form  the  roof  shield  three  slabs  were  laid  side  by  side  so  that 
the  stepped  sides  overlapped,  thus  eliminating  the  possibility  of 
radiation  streaming  between  the  slabs.  A  thin  layer  of  concrete  was 
placed  on  the  top  surface  of  the  silo  wall  immediately  prior  to  the 
positioning  of  the  first  layer  of  roof  slabs.  This  ensured  a  complete 
seal  between  the  roof  and  silo  wall  thereby  eliminating  any  possibility 
of  streaming  under  the  lowest  layer  of  slabs.   To  minimize  the  effect 
of  the  steel  frames  on  the  penetration  of  radiation  through  the  roof, 
the  slabs  were  positioned  such  that  the  long  axis  of  the  slabs  were 
directed  toward  the  measurement  locations.  With  this  configuration  the 
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frames  would  significantly  affect  only  that  portion  of  the  field 
leaving  the  silo  perpendicular  to  the  source/detector  axis. 

The  thickness  of  each  slab  was  measured  at  six  locations  with  a 
pair  of  specially  constructed  long  armed  capilers.   The  results  are 
shown  in  Table  1  .   In  order  to  provide  a  consistent  geometry,  the 
slabs  were  always  positioned  in  the  same  stacking  pattern  as  shown  in 
Fig.  15. 

Table  1.  Mean  Thicknesses  of  Silo  Roof  Shields** 


Slab 

No. 

Thickne 

ss  (cm) 

1 

10.77 

± 

0.29 

2 

10.90 

± 

0.48 

3 

10.72 

± 

0.19 

4 

10.40 

+ 

0.16 

5 

10.34 

+ 

0.21 

6 

11.03 

+ 

0.34 

7 

10.45 

± 

0.25 

8 

10.53 

+ 

0.19 

9 

10.45 

+ 

0.16 

10 

10.40 

± 

0.22 

11 

10.58 

+ 

0.16 

12 

10.58 

± 

0.24 

Center,  9 

+ 

12 

21.03 

+ 

0.29 

Center,  1+6 

+ 

9  +  12 

42.83 

+ 

0.53 

Standard  concrete  test  specimens  were  taken  for  ail  of  the  con- 
crete pours  of  the  silo  wall,  collimator  wedges,  and  roof  slabs.   The 
densities  of  these  specimens  were  determined  by  faculty  members  of  the 
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Civil  Engineering  Department  at  Kansas  State  University.   These  results 

Q 

are  tabulated  in  Table  2  ,  where  no  correction  has  been  made  for  the 
affect  of  the  structural  steel  within  the  actual  components. 


Table  2.  Measured  Concrete  Densities  for  Various  Source 
Silo  Structures8 

3 
Structure  Density  (g/cm  ) 


Silo  2.17  ±  0.02 

Roof  Shields  2.13  ±  0.02 

Collimator  Wedges  1.93  ±  0.02 


2.3.4  Source  Irradiators 

The  source  irradiators  were  designed  and  manufactured  by  J.  L. 
Shephard  and  Associates,  Inc.   These  irradiators  were  mounted  on  a 
moveable  platform  in  such  a  way  that  either  cask  could  be  positioned 
at  the  center  of  the  silo.   Plates  9  and  10  shows  the  two  irradiator 
units  during  their  installation  while  Fig.  16  gives  the  dimensions  of 
the  irradiator  assembly.   The  larger  cask,  contains  an  effective  3S00  Ci 
Co-60  source  while  the  smaller  cask  contains  both  the  effective  10.3 
and  229  Ci  Co-60  sources*.   The  two  largest  sources  when  in  the  exposed 
position  had  their  midpoint  the  apex  of  the  collimation  cone  while  the 
10.3  Ci  source  was  2"  from  the.  aoex. 


The  source  strengths  as  reported  by  the  manufacturer  were  12.2,  255,  and 
5220  Ci.   To  ascertain  the  effective  source  strengths,  a  "good  geometry" 
experiment  was  performed  at  a  distance  of  300  m.   These  results  were  then 
used  to  cross  calibrate  the  skyshine  sources  using  a  N3S  source  of  known 
intensity  as  a  reference.   The  observed  reduction  in  intensity  is  pre- 
sumbably  due  to  shelf-shielding  effects. 


4ft 


Plate  9.   Source  irradiator  prior  to  placement  in  source  silo. 
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The  irradiators  can  be  operated  from  a  remote  console  up  to  100  m 
from  the  sources.  During  operation,  the  central  rectangular  drawer, 
or  attenuator,  at  the  top  of  each  cask  is  withdrawn  pneumatically  until 
a  vertical  hole  in  the  slide  is  properly  positioned  over  the  appropriate 
vertical  source  hole  in  each  cask.   Of  the  two  source  holes  in  the 
smaller  cask,  the  one  for  the  229  Ci  source  is  on  the  center  line  of 
the  unit  while  the  10  Ci  source  position  is  2"  from  that  point.  After 
the  attenuator  is  accurately  positioned,  the  specified  source  is  raised 
pneumatically  out  of  the  irradiator  to  the  "exposed"  position.  All 
movements  of  the  irradiator  systems  are  monitored  by  microswitches  and 
status  lights  on  the  operator's  console  indicate  the  position  of  the 
sources. 

To  protect  the  entire  irradiator  assembly,  a  sheet  metal  weather 
housing  was  placed  over  both  units  and  secured.   Incorporated  into  the 
metal  cover  were  two  thin  walled  hemispherical  steel  domes  positioned 
over  the  two  source  locations.   Therefore  regardless  of  the  upward  direction 
taken  by  the  gamma  photons,  all  radiation  would  experience  identical 
amounts  of  mass  attenuation  while  penetrating  the  weather  cover. 

2.4  Measurement  Procedures 

2.4.1  Trailer  Positioning 

For  each  skyshine  spectral  measurement,  it  was  required  that  the 
detector  cavity  axis  be  horizontal  and  directed  toward  the  silo.   In 
order  to  adjust  the  attitude  of  the  trailer,  an  I-beam  frame  was  welded 
beneath  the  rear  of  the  trailer  so  that  hydraulic  jacks  could  be  placed 
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under  each  rear  corner.   In  this  way  the  entire  back  of  the  trailer  could 
be  raised  and/or  tilted.  Also,  by  simply  adjusting  the  length  of  the 
front  support  legs,  the  front  end  could  be  raised  or  lowered. 

During  the  construction  of  the  collimator  assembly,  the  trailer 
was  positioned  using  the  above  adjustment  procedures  so  that  the  lower 
edge  of  the  shadow  shield  was  horizontal  and  the  sides  vertical.  Two 
sets  of  brackets  were  then  securely  fastened  to  the  rear  and  one  side 
of  the  trailer  so  that  when  a  level  rested  on  them  a  level  condition  was 
indicated.  Thus  for  routine  positioning,  adjusting  the  attitude  of  the 
trailer  until  the  two  perpendicular  lines,  defined  by  the  two  sets  of 
brackets,  defined  a  horizontal  plane  ensured  that  the  collimator  was 
vertical  and  the  detector  assembly  axis  horizontal. 

In  order  to  align  the  detector  assembly  with  the  source  position, 
a  transit  was  used  to  place  two  aiming  pins  on  opposite  ends  of  the 
peak  of  the  environmental  shed  roof  so  that  they  defined  a  vertical 
plane  with  the  cavity  axis.   Therefore,  by  positioning  the  trailer 
so  that  the  two  aiming  pins  and  the  center  of  the  source  silo  were 
colinear  ensured  that  the  detector  cavity  axis  was  also  aligned  with 
the  source  silo.   Since  by  necessity  this  alignment  was  done  before 
the  trailer  was  leveled,  some  error  in  alignment  might  have  occurred 
although  it  was  estimated  to  be  less  than  2°.   It  was  demonstrated  by 
experiment  that  the  measured  skyshine  spectrum  was  totally  insensitive 
to  such  small  errors  in  the  trailer  orientation. 
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2.4.2  Spectral  Quality  Control 

To  ensure  that  the  detection  system  was  properly  calibrated, 
several  test  spectra  were  taken  prior  to  any  skyshine  measurements. 
After  the  instrumentation  had  been  turned  on  for  approximately  30  minutes, 
a  mixed  isotope  source  was  used  to  adjust  the  system  to  the  proper  cali- 
bration.  Next  a  Cs-137  response  function  was  taken  which  was  analyzed 
with  respect  to  the  total  count  rate  as  a  check  on  the  system  efficiency. 
This  spectrum  was  also  inspected  visually  to  verify  the  spectral  shape. 
Once  it  was  confirmed  that  the  system  was  operating  properly,  a  back- 
ground spectrum  was  taken  after  which  skyshine  measurements  were  made. 
After  each  skyshine  spectrum,  a  mixed  isotope  response  function  and 
a  background  spectrum  were  taken  in  order  to  monitor  the  system  cali- 
bration and  background  levels  with  time. 

The  count  times  for  each  skyshine  spectral  measurement  were  chosen 
so  that  acceptable  statistical  uncertainties  in  the  measured  spectra 
would  be  obtained.   Additionally,  duplicate  measurements  were  made  to 
determine  the  reproducibility  of  the  results  and,  when  possible,  different 
sources  were  used  as  further  test  of  reproducibility. 
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3.0  Spectral  Analysis  and  Comparisons 

A  major  aspect  of  the  Benchmark  Skyshine  Experiment  was  the  measure- 
ment of  skyshine  energy  spectra  with  a  collimated  Nal(Tl)  spectrometer. 
One  of  the  reasons  for  the  measurement  of  the  skyshine  spectra  was  to 
serve  as  a  benchmark  against  which  the  results  of  photon  transport  models 
could  be  compared.  However,  because  of  the  statistical  nature  of  the 
energy  deposition  mechanisms  in  the  crystal  and  the  subsequent  pulse 
shaping  and  amplification  processes,  the  observed  pulse-height  spectrum 
often  bears  little  semblance  to  the  true  incident  spectrum.   It  is  thus 
necessary  to  correct  or  "unfold"  the  measured  spectra  to  obtain  the  actual 
energy  distribution  of  the  incident  radiation.  The  unfolding  technique 
used  in  this  study  to  obtain  the  energy  dependence  of  the  original  spectrum 
from  the  observed  spectrum  will  be  described  in  the  first  section  of  this 
chapter.   The  second  section  of  this  chapter  compares  the  unfolded  measured 
spectra  with  the  energy  spectra  calculated  with  the  DOT-3.5"1  photon  trans- 
port code. 

3.1  Data  Analysis  Techniques 

A  common  problem  in  an  experimental  measurement  is  the  distortion 
of  the  quantity  of  interest  by  the  detection  device.   In  the  case  of 
measuring  the  energy  distribution  of  a  radiation  field  with  a  Nal  de- 
tector, the  differences  between  the  incident  and  measured  spectrum  arise 
in  the  crystal  as  well  as  in  the  electronics  of  the  system.  The  crystal 
induced  distortions  arise  mainly  from  the  statistical  nature  of  the 
gamma-ray  interactions  within  the  crystal  resulting  in  either  total  or 
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partial  deposition  of  the  incident  photon  energy.   The  primary 
component  of  the  electronic  distortion  is  the  statistical  fluctuation 
of  the  photoelectron  multiplication  in  the  photomultiplier  tube  (PMT) .  To 
obtain  an  estimate  of  the  true  energy  dependence  of  the  primary  radiation 
field,  these  distortions  must  be  corrected  by  "unfolding"  the  measured 
spectrum. 

3.1.1  Spectral  Characteristics 

The  height  of  the  voltage  pulse  entering  the  multichannel  analyzer 
(MCA)  is  directly  proportional  to  the  amount  of  energy  deposited  in  the 
crystal  by  the  incident  photon.  Although  several  types  of  interactions 
are  possible  between  the  photon  and  the  crystal,  three  types  play  the 
dominant  role  in  radiation  spectroscopy:   photoelectric  absorption, 
Compton  scattering,  and  pair  production  which  predominate  below  .3  MeV, 
between  .3  and  1.1  MeV,  and  above  1.1  MeV  respectively  .  All  of  these 
processes  result  in  either  the  partial  or  complete  transfer  of  the 
photon  energy  to  a  secondary  electron. 

An  understanding  of  the  energy  deposition  and  signal  collection, 
shaping  and  amplification  processes  in  a  Nal  spectrometer  form  a  basis 
for  interpreting  the  observed  pulse  height  distribution.   Of  special  con- 
cern are  the  following  processes  which  introduce  non-ideal  responses  in 

9 
the  spectrometer  system  : 

1)  partial  energy  deposition  of  the  incident  photon  energy  in 
the  crystal, 

2)  nonlinear  light  production  within  the  crystal, 

3)  non-uniformity  of  photoelectron  collection  from  the  photo- 
cathode,  and 
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4)   statistical  fluctuations  in  the  electron  multiplication  during 
the  cascade  through  the  dynodes. 

A  typical  Cs-137  response  function  is  shown  in  Fig.  17,  where 
some  of  the  more  obvious  features  have  been  identified.   The  full  energy 
peak,  or  photopeak,  arises  when  the  entire  energy  of  the  incident 
photon  is  deposited  in  the  crystal  within  a  time  interval  too  short 
to  be  resolved  by  the  detector.   The  deposition  may  occur  in  one  photo- 
electric absorption  or  a  series  of  Compton  scatters  followed  by  a 
photoelectric  absorption. 

One  should  notice  that  even  though  the  photopeak  represents  the 
detection  of  a  photon  with  a  discrete  energy,  a  Gaussian  shaped  response 
is  obtained.   This  is  the  result  of  many  factors,  the  dominant  one  being 
the  statistical  nature  of  the  electron  multiplication  cascade  process  in 
the  PMT.  The  multiplication  factor  for  a  given  dynode,  defined  as  the 
number  of  secondary  electrons  emitted  per  incident  electron,  is  not  con- 
stant but  fluctuates  about  some  mean  value.   This  in  turn  causes  the 
amplitude  of  the  PMT  output  voltage  pulse  (produced  by  the  deposition  of 
energy  E  in  the  crystal)  to  fluctuate  about  some  mean  value.   As  a  conse- 
quence of  this,  the  photopeak  is  distributed  about  some  mean  channel 
corresponding  to  energy  E.   This  electronic  smearing  is  obviously  not 
limited  to  the  photopeak  but  is  manifested  over  the  entire  spectrum. 

Another  contributing  factor  to  the  spectral  broadening  is  non- 
linear crystal  response  with  energy.   If  the  light  production  within 
the  crystal  varies  with  the  energy  of  the  electron,  the  total  light 
yield  will  be  different  from  event  to  event,  depending  on  how  the 
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energy  of  the  incident  photon  is  subdivided  between  the  various  secondary 
electrons.  This  effect  would  be  expected  to  be  significant  for  initial 

Q 

energies  above  400  keV  where  multiple  interactions  are  predominant  . 
Some  additional  broadening  occurs  due  to  nonuniform  photoelectron 
production  at  the  photocathode  and  the  subsequent  collection  at  the 
first  dynode. 

The  broadening  of  the  lower  side  of  the  photopeak  is  influenced 
by  the  escape  of  the  secondary  photon  after  undergoing  two  or  more 
Compton  interactions.   Shallow  angle  scattering  of  the  incident  photons 
before  reaching  the  detector  is  another  contributing  factor.   The  shape 
of  the  high  energy  side  of  the  photopeak  is  affected  in  high  count  rate 
situations  by  the  superposition  of  one  or  more  pulses  at  the  photocathode. 
This  distortion  is  commonly  referred  to  as  "pile-up". 

The  Compton  distribution,  or  plateau,  represents  a  series  of 
events  in  which  scattered  photons  escape  the  crystal  without  depositing 
all  of  their  initial  energy  in  the  detector.   Since  the  energy  of  the 
secondary  photon  is  a  continuous  function  of  the  scattering  angle,  the 
Compton  plateau  extends  from  zero  energy  up  to  the  energy  corresponding 
to  the  maximum  amount  of  energy  that  can  be  transferred  to  an  electron 
in  a  single  event.   This  cutoff  is  commonly  called  the  Compton  edge  and 
corresponds  to  a  180°  angle  scattering  event.   In  practice,  the  Compton 
edge  is  not  a  sharp  cutoff  because  of  shallow  angle  and  multiple 
scatters  which  tend  to  fill  in  the  "valley"  between  the  Compton  edge  and 
the  photopeak. 
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The  backscatter  peak  is  caused  by  those  photons  which  have  ex- 
perienced a  large  angle  scattering  event  prior  to  reaching  the  crystal. 
As  can  be  seen  in  Fig.  18,  the  energy  of  the  recoil  photon  changes  very 
slightly  for  scattering  angles  greater  than  120°.   Therefore,  a  peak  is 
formed  as  a  result  of  those  photons  depositing  all  of  their  energy  in 
the  crystal  and  is  centered  around  the  channel  corresponding  to  the 
scattered  photon  energy  after  a  180°  Compton  scatter. 

The  characteristic  x  ray  also  shown  in  Fig.  17,  is  a  result  of  an 
internal  conversion  process  within  the  daughter  of  the  isotope  under 
consideration.  As  such,  this  process  in  itself  represents  a  monoenergetic 
photon  source  and  although  sometimes  proving  useful  in  identification, 
often  complicates  the  interpretation  of  the  response  of  the  system  to 
the  primary  gamma  emission  at  very  low  energies. 

A  very  low  energy  "iodine  escape  peak"  is  often  noticed  in  a  Nal 
pulse  height  distribution  in  association  with  low  energy  photon  sources. 
During  the  energy  deposition  process,  iodine  atoms  within  the  Nal  crystal 
are  excited  and  emit  x  rays  at  approximately  29  keV  upon  reversion  to  the 
ground  state.   When  dealing  with  low  energy  sources  such  as  an  internal 
conversion  x  ray,  the  point  of  interaction  is  usually  close  to  the  sur- 
face of  the  crystal  since  the  attenuation  coefficient  for  Nal  increases 
significantly  at  low  energies.   Being  produced  near  the  periphery  of  the 
crystal,  the  de-excitation  x  rays  have  a  high  probability  of  escaping 
the  crystal.   The  escape  of  the  x  ray  causes  the  energy  deposited  in  the 
crystal  to  be  29  KeV  below  the  energy  of  the  incident  photon  and  an 
iodine  escape  peak  is  formed. 
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Scattering  Angle 
Fig.  18.   Secondary  photon  energy  as  a  function  of 
scattering  angle. 
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Two  other  escape  peaks  may  be  encountered  when  dealing  with  source 
energies  above  1.1  MeV.   As  a  consequence  of  pair  production,  two  0.51 
MeV  annihilation  photons  are  produced.   One  or  both  of  these  photons 
may  escape  the  crystal  without  an  interaction  resulting  in  a  first-  or 
second-escape  peak,  respectively. 

3.1.2  Theory  of  Unfolding 

The  pulse  height  distribution  of  the  output  voltage  pulses  of  a 
scintillation  spectrometer,  C(V),  is  related  to  the  incident  gamma-ray 

fluence,  <t(E'),  by 

E 
max 

dE'  R(E',V)*(E')  =  C(V)  (3.1-1) 

0 

where  R(E' ,V)  is  the  probability  that  an  incident  photon  of  energy 

E'  will  produce  an  output  pulse  within  unit  voltage  about  V,  and  E 

max 

is  the  maximum  energy  of  the  primary  spectrum.   The  function  R  is  often 
referred  to  as  the  response  surface  of  the  spectrometer  because  it  describes 
the  response  of  the  system  to  photons  of  a  given  incident  energy,  E'.  The 

binning  or  discretization  process  performed  by  the  MCA  can  be  written  as 

E 
(■     f  max 

C.  -     dV     dE'  R(E\V)*(E')    1-1,2,.  -N  (3.1-2) 

AV.     ° 

E 
I   max 

C   =     dE'  *(E')R ,(E')  1-1,2,.. .8  (3.1-3) 

'o  X 


where 


f 

ME')  -  J   dVR(E'.V)  1=1,2, ...N,         (3.1-4) 

AV. 
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is  called  the  "spectrometer  response  function"  and  physically  may 

be  interpreted  as  the  expected  number  of  counts  observed  in  channel  i 

of  the  pulse  height  distribution  due  to  a  unit  fluence  of  photons  of 

energy  E1.   In  this  formulation  N  is  the  number  of  channels  in  the  pulse 

height  spectrum  and  C  is  the  observed  number  of  counts  in  channel  i 

corresponding  to  the  observed  number  of  voltage  pulses  ranging  from 

V.  to  V.+AV.  in  amplitude. 
1     i   1      r 

The  response  surface  of  the  spectrometer,  a  typical  representation 
of  which  is  shown  in  Fig.  19  ,  can  be  determined  in  a  variety  of  ways. 
Much  effort  has  been  expended  in  using  Monte  Carlo  techniques  to  pre- 
dict the  response  of  a  specific  detection  system  to  a  source  of 

.        t     11-13 
monoenergetic  photons     .  Although  excellent  results  have  been 

obtained,  the  system  response  is  highly  sensitive  to  parameters  such 

as  detector  size  and  collimator  geometry  which  necessitates  independent 

calculations  for  each  detection  system.   This  procedure  is  prohibitive 

in  most  cases  because  of  the  large  computational  effort  involved. 

A  second  technique  involves  experimentally  measuring  the  response 
of  the  system  at  several  discrete  energies  and  using  elaborate  inter- 
polation schemes  to  calculate  the  response  at  those  energies  which  are 
not  available  experimentally14'15.   This  technique  is  widely  used  al- 
though the  source  preparation,  measurements,  background  corrections  and 
calibration  matching  also  make  this  technique  quite  onerous  to  use. 

To  obtain  the  actual  energy  distribution  of  the  incident  photon 
fluence  from  the  measured  pulse-height  distribution,  it  is  convenient 
to  restate  Eq.  (3.1-3)  in  matrix  notation  as 
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MONOENERGETIC     GAMMA- BAY     SOUBCES 
USED     TO    DETERMINE    RESPONSE 


Fig.  19....  Typical  spectrometer  response  surface.   (From 
R.  L.  Heath,  Ref.  10) 


62 


£"&*  (3.1-5) 

where  C  is  the  number  of  counts  in  channel  i,  R.  .  represents  the 
i  XJ 

number  of  counts  in  channel  i  per  unit  fluence  in  energy  interval  j, 

and  $     represents  the  incident  photon  fluence  in  energy  interval  j . 

Although  direct  solution  of  Eq.  (3.1-5)  for  *  is  possible 

J 

using  standard  inversion  techniques,  unphysical  oscillations  and 

14 
negative  values  often  appear  in  the  results   .   Therefore,  various 

iterative  schemes  have  been  devised  in  an  effort  to  consistently  obtain 

a  stable,  physically  significant  solution     . 

3.1.3  Derivation  of  the  SEGO  Method 

The  technique  used  in  this  study  to  unfold  the  measured  skyshine 
spectra  was  developed  by  Young  and  Burrus  and  is  designed  to  correct 
for  the  distortions  in  the  spectrum  due  to  partial  energy  deposition  in 
the  crystal  while  at  the  same  time  preserving  the  inherent  spectral 
smoothing  caused  by  the  finite  resolution  of  the  detector  system.   In 
this  technique,  the  Compton  plateau  and  the  various  scattering  and 
escape  peaks  are  "stripped"  from  the  observed  pulse-height  distribution. 
The  resulting  spectrum  is  an  approximation  to  the  incident  radiation 
field  which  has  been  broadened  according  to  the  system  resolution. 

The  observed  pulse  height  distribution  is  the  result  of  two 
distinct  processes,  the  first  being  the  deposition  of  energy  within 
the  crystal,  the  second  being  the  broadening  effect  in  converting  the 
deposited  energy  into  a  voltage  pulse.   The  energy  deposition  process 
can  be  formulated  in  terms  of  a  function  U  defined  as 
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U(E,E')dE' 


the  probability  that  an  incident  photon  of 
energy  E  will  deposit  energy  in  the  crystal 
within  dE'  about  E' 


Similarly,  the  spectral  broadening  can  be  expressed  by 

(the  probability  that  the  deposition  of  energy  E  in 
S(E,E')dE'  =  the  crystal  will  produce  a  PMT  voltage  pulse 

(corresponding  to  a  photon  with  energy  in  dE'  about  E1 

Physically,  the  incident  spectrum  due  to  a  monoenergetic  source 

is  ideally  represented  by  a  delta  function  at  some  source  energy  E.  This 

energy  is  deposited,  either  in  part  or  in  full,  in  the  crystal.   The 

response  of  the  crystal  is  then  broadened  by  the  trailing  electronics. 

The  assumption  made  in  the  development  of  the  SEGO  method  is  that  the 

pulse  height  distribution  obtained  from  the  above  sequence  of  events  is 

equivalent  to  what  would  be  observed  if  the  incident  spectrum  was  first 

smeared  by  the  inherent  resolution  of  the  spectrometer  and  then  detected 

by  a  detection  system  with  perfect  resolution.  With  this  alternative 

viewpoint  the  spectrometer  reponse  function,  R. (E'),  of  Eq.  (3.1-4) 

can  be  expressed  in  terms  of  U  and  S  as 

E 
t  f   max 

St(E')  =     dE     dE"  U(E",E)  S(E',E")   1=1, 2,... N  .      (3.1-6) 
AE.     0 

It  should  be  pointed  out  that  in  this  formulation  the  smearing  function 

S,  is  a  function  of  the  source  energy  E'  and  not  of  the  energy  deposited 

in  the  crystal,  E.   Equation  (3.1-6)  can  be  rewritten  as 
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R^X')  -  j  dE  I       J  dE"  U(E",E)  S(E',E")   i=l,2,...N.      (3.1-7) 
AE.    k=1  ML 

If  one  assumes  that  U  and  S  vary  slowly  across  any  energy  bin  AIL  or 

AE^  Eq.  (3.1-7)  becomes 

N 
R.(E')  =  I   U(Ek,Ei)AEkAE±  S(E'  ,Bj,)   1=1,2,. ..N  .  (3.1-8) 

With  the  definition 

Uik  =  AEk  AEi  U(Ek'V  '  "•1"9> 

Eq.  (3.1-8)  can  be  expressed  as 

R.(E')  =  I   U±k  8(8', E.)   i=l,2,...N  .  (3.1-10) 

k 

Substitution  of  Eq.  (3.1-10)  into  Eq.  (3.1-3)  yields 

N      i   max 
C  =  I   U      dE'  *(E')  S(E',IL)    i-1,2, ...N  .  (3.1-11) 

k-1     '0  K 

Notice  that  the  quantity  under  the  integral  is  merely  the  incident 
fluence  which  has  been  smeared  according  to  the  system  resolution. 

One  can  therefore  define 

E 
t   max 

<S(Ek)  =  j  dE'  $(E')  S(E',Ek)  (3.1-12) 

as  the  broadened,  incident  spectrum.   Substitution  of  this  definition 
into  Eq.  (3.1-11)  then  yields 

°i  =  I   Uik  *(Ek)     i=l,2,...N  .  (3.1-13) 
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It  is  the  discretized,  smeared  incident  spectrum,  *(E, ),  that  is 
solved  for  in  the  SEGO  code. 

Equation  (3.1-13)  represents  a  set  of  N  equations  in  N  unknowns, 
namely  the  $(K)'s.      The  matrix  elmenets  U  ,  are  a  discrete  repre- 
sentation of  the  response  surface  of  the  ideal  spectrometer  (i.e.,  a 
spectrometer  with  perfect  resolution).   Physically,  U.,  may  be 
interpreted  as  the  probability  a  count  will  be  recorded  in  channel  i 
per  incident  photon  on  the  Nal  crystal  whose  energy  corresponds  to 
channel  k.   The  realistic  assumption  is  made  that  it  is  impossible  for 
a  photon  of  energy  E"  to  deposit  an  energy  E  in  the  crystal  with  E  >  E". 
Therefore,  U   is  identically  zero  for  those  values  of  i  >  k.  However, 
this  assumption  ignores  the  effects  of  sum  peaks  where  two  photons  deposit 
their  entire  energy  in  a  time  interval  too  short  to  be  resolved  by  the 
crystal.   Also  ignored  are  pile-up  effects. 

Under  this  assumption,  the  system  of  equations  defined  by  Eq. 
(3.1-13)  becomes  triangular  in  nature.   Denoting  $(K  )  by  X  ,  Eq.  (3.1-13) 
can  be  written  as 


11  1 


U12X2+U13X3+-- 

••  +Disxs-Ci 

u12x2  +  u23x3  +  .. 

••    +U2NXN-C2 

u33x3  +  .. 

-    +U3NXN=?3 

umh  =  V 

(3.1-14) 


This  system  can  be  immediately  solved  by  first  solving  for  X 
and  working  backwards  to  X  .   In  general,  the  solution  can  be  written  as 
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Si  -  K  "  J+1D*n  XJ/Utr  1=N"1.N-2 1  •    (3.1-15) 

This  solution  (which  is  related  to  the  true  incident  spectrum  by 
Eq.  (3.1-12))  gives  an  approximate  representation  of  the  incident 
fluence  on  the  Nal  crystal  which  is  corrected  for  those  photon-crystal 
interactions  which  do  not  deposit  all  of  the  photon's  energy  in  the 
crystal,  i.e.,  only  interactions  leading  to  counts  in  the  photopeaks 
are  retained.  Moreover,  the  system  resolution  has  been  imposed  on  this 
unfolded  spectrum  since  it  is  fruitless  to  try  to  resolve  energy  structure 
outside  of  the  inherent  resolution  of  the  spectrometer. 

In  the  preceding  unfolding  description,  the  statistical  uncertainty 
in  the  observed  values  of  C.  has  been  ignored.   However,  due  to  the 
statistical  nature  of  the  detection  process,  the  measured  distribution 
is  only  one  member  of  a  family  of  possible  results.   A  stochastic 
perturbation  technique  is  used  in  the  SEGO  code  to  evaluate  the 
uncertainty 'in  the  unfolded  result.   This  method  eliminates  the  need 
for  the  inversion  of  the  entire  response  matrix  which  would  be  required 
if  standard  error  propagation  techniques  were  used.   Given  a  measured 

pulse-height  distribution,  C  ,  i=l,2, N,  each  with  a  corresponding 

standard  deviation,  S.,  i=l,2,...N,  one  can  stochastically  perturb  the 
given  spectrum  by  adding  to  each  C.  a  randomly  selected  normal  deviate 
chosen  from  the  Gaussian  distribution  with  mean  zero  and  standard  deviation 
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4 
S  .   By  generating  and  unfolding  several  of  these  perturbed  spectra, 

the  standard  deviation  of  the  unfolded  measured  spectrum  can  be  esti- 
mated from  the  range  of  the  results.   In  the  SEGO  code,  three  spectra 
generated  by  perturbing  each  channel  in  the  above  fashion,  are  unfolded 
along  with  the  original  observed  spectrum. 

3.1.4  Description  and  Construction  of  the  Response  Matrix 
The  distinguishing  characteristic  of  the  SEGO  method  is  that  each 
response  function,  U    is  constructed  through  the  superposition  of 
appropriately  placed  rectangles  as  suggested  by  Hyodo  .   This  concept 
is  illustrated  in  Fig.  20,  where  one  column  of  the  U  matrix  is  shown 
in  the  idealized  form.  Notice  that  each  feature  of  the  response  function 
is  approximated  by  a  discrete,  sharp  component.  A  long  rectangle  ex- 
tending from  zero  energy  to  the  Compton  edge  is  used  to  approximate  the 
Compton  plateau.   The  "valley  fill  function",  located  immediately  below 
the  photopeak  is  used  to  approximate  those  counts  which  occur  between 
the  full  energy  peak  and  the  Compton  edge.  Other  features  of  the 
response  function  include  the  representation  of  the  first-  and  second- 
escape  peaks  and  the  backscatter  peak.   Additional  features  could  be 
added  to  account  for  other  minor  structures  if  necessary. 

As  is  obvious  from  Fig.  20,  large  portions  of  the  idealized  response 
functions  are  flat.   This  flatness  is  used  in  the  SEGO  code  to  reduce 
significantly  the  execution  time.   Returning  to  Eq.  (3.1-14)  and  sub- 
tracting equation  (j+1)  from  the  j-th  equation  for  j=l...N,  one  obtains 
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*uh  +  VwV  +  h^u-itf  +  •••  +  Vuin"V  -  crc2 

U22X2       +  X3(U23-U33)  +  ...  +  ^(B^-V  ■  C2-C3 

u33x3   +...+x^(u3N-u4N)  =  c3-c4 

With  this  transformation  of  the  system,  many  coefficients  are  now 

zero.  The  only  nonzero  components  are  located  where  the  U,,  's  have 

ik 

discontinuities.  Thus  the  transformed  U  matrix  is  still  triangular 
and  now  is  only  sparsely  filled.   The  summation  in  Eq.  (3.1-15)  needs 
only  to  be  carried  out  over  the  nonzero  components  of  which  there  are 
only  ten:   the  upper  and  lower  edges  of  the  backscatter  peak,  the  first- 
and  second-escape  peaks  and  the  "valley  fill"  function,  plus  the  Compton 
edge  and  the  photopeak.   This  transformation  is  reversed  during  the 
solution  of  the  system. 

One  of  the  advantages  of  the  Hyodo  response  functions  is  that 
the  spectrum  features  vary  fairly  slowly  with  energy.   Thus  it  is 
feasible  to  tabulate  them  at  a  few  energies  and  use  linear  interpolation 
to  obtain  the  necessary  values  for  other  energies  during  the  course  of 
the  computation.   This  property  eliminates  the  need  to  store  the  entire 
response  matrix. 

Typically,  eight  parameters  are  needed  to  define  a  Hyodo-type 
response  function  for  an  incident  monoenergetic  photon  source  of  energy 
E' .  They  are: 
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1)  EF,  defined  as  the  ratio  of  the  total  counts  in  the  spectrum 
to  the  number  of  photons  reaching  the  crystal  face, 

2)  PT,  defined  as  the  ratio  of  the  area  of  the  photopeak  to  the 
area  of  the  total  spectrum  (peak  to  total  ratio) , 

3)  ES,  defined  as  the  first-escape  peak  to  total  ratio, 

4)  EC,  defined  as  the  second-escape  peak  to  total  ratio, 

5)  BK,  defined  as  the  backscattering  peak  to  total  ratio, 

6)  VH,  defined  as  the  valley  to  total  ratio, 

7)  TO,  defined  as  the  valley  upper  edge  energy,  and 

8)  VL,  defined  as  the  valley  width  in  energy  units. 

Experimental  response  functions  were  measured  using  Am-241, 
Hg-203,  Cs-137,  and  Co-60  calibration  sources  which  were  positioned 
on  the  detector  axis  at  a  distance  of  69.85  cm.  Three  response 
functions  from  each  source  were  measured  to  ensure  reproducibility. 
Long  count  times  were  also  used  in  order  to  minimize  the  statistical 
uncertainties  in  the  result.  An  appropriate  background  spectrum  was 
subtracted  from  each  experimental  response  function  before  analysis. 

An  initial  response  matrix  was  calculated  from  these  70-cm  experi- 
mental response  functions.   Since  all  of  the  sources  emitted  photons 
of  two  different  energies,  some  approximations  were  necessary  in 
order  to  treat  them  as  monoenergetic  sources.   For  example,  in  the 
case  of  Cs-137,  the  area  of  the  32.88  keV  x  ray  must  not  be  included 
in  the  calculation  of  the  unfolding  parameters  for  the  .662  MeV  photon. 
To  this  end,  the  x  ray  peak  was  eliminated  by  assuming  that  the  Compton 
plateau  was  continuous  through  the  region  beneath  the  peak.   Since  there 
should  be  very  little,  if  any,  contribution  to  the  Compton  continuum 
from  the  x  ray,  the  resulting  spectrum  should  very  closely  approximate 
the  response  of  the  system  to  a  monoenergetic  source  of  0.662  MeV  photons. 
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In  the  case  of  Co-60,  the  approximation  described  was  not  valid 
since  source  photons  of  both  energies  (1.173  and  1.332  MeV)  make  a 
significant  contribution  to  Compton  plateau.   Therefore,  the  experi- 
mental Co-60  response  function  was  treated  as  monoenergetic  source 
with  an  energy  of  1.253  MeV. 

The  valley  upper  edge  was  determined  from  the  lower  side  of  the 
photopeak.   The  sharply  increasing  region  of  the  peak  was  used  to 
define  a  line  which  would  intersect  at  the  approximate  lower  edge  of 
the  photopeak.  The  valley  lower  edge  was  at  first  taken  to  be  the 
Compton  edge,  but  was  later  adjusted  to  extend  slightly  below  this 
point. 

The  initial  results  obtained  from  unfolding  the  four  calibration 
spectra  with  the  70-cm  parameters  were  not  acceptable  so  individual 
values  were  adjusted  by  trial  and  error  until  most  of  the  non-photopeak 
features  were  eliminated  and  the  photopeak  area  was  within  a  few  percent 
of  the  expected  value.   The  unfolded  70-cm  response  functions  are  shown 
in  Fig.  21,  where  although  the  various  secondary  structures  have  not 
been  completely  eliminated,  the  area  outside  of  the  photopeak  is  less 
than  V/.   in  all  cases.   The  optimized  70-cm  unfolding  parameters  are 
shown  in  Table  3. 

Because  the  effective  source/detector  distance  for  the  skyshine 
incident  flux  would  be  much  greater  than  70-cm,  a  significant  change 
in  the  total  or  intrinsic  efficiency  of  the  detector  would  be  experienced. 
The  intrinsic  efficiency,  e  (S),  defined  as  the  probability  that  a 
photon  incident  upon  the  detector  face  will  interact  within  its  active 
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Table  3.      Optimized   70-cm  Unfolding  Parameters 
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E (MeV) 

EF 

PT 

ES 

EC 

BK 

VH 

VU 

VL 

0.0 

1.113 

0.9525 

0.0 

0.0 

0.0070 

0.0 

0.0 

0.0 

0.2 

1.083 

0.8304 

0.0 

0.0 

0.0110 

0.1069 

0.1738 

0.0859 

0.4 

1.052 

0.7239 

0.0 

0.0 

0.0160 

0.0958 

0.3610 

0.1240 

0.6 

1.0216 

0.6312 

0.0 

0.0 

0.0210 

0.0898 

0.5481 

0.1620 

0.8 

.9912 

0.5502 

0.0 

0.0 

0.0250 

0.0858 

0.7352 

0.2001 

1.0 

.9607 

0.4797 

0.0 

0.0 

0.0300 

0.0828 

0.9224 

0.2381 

1.2 

.9301 

0.4182 

0.005 

0.004 

0.0340 

U.0804 

1.1095 

0.2762 

1.4 

.8996 

0.3646 

0.010 

0.008 

0.0390 

0.0785 

1.2967 

0.3142 

1.6 

.8691 

0.3179 

0.015 

0.012 

0.0440 

0.0768 

1.4838 

0.3523 

75 

19 
volume,  can  be  calculated  as 

J   (l-e"UX)  dn 

Q 

et<B)  -  °    g                             (3 

.1-17) 

0 

where  fi  is  the  solid  angle  subtended  by  the  crystal. 

In  the  case  of  a  point  source  located  on  the  axis  of  a  right 

cylindrical  detector  of  height  t  and  radius  r,  as  illustrated  in 

Fig.  22,  Eq.  (3.1-17)  becomes 

e 

fm   sineu-e-^^) 

«eW  -  s°                            (3 

1-18) 

J  d6  sin9 

0 

where 

•~S                for  0  <  8  <   9  ' 
cosS             —   — 

x(8)  = 

(3 
r-h-tan9       ., 

1-19) 

.  n     tor  o  <  d  <  o   , 
sino          —   —   m 

8  =  tan   (r/h),                                  (3 

1-20) 

9'  ■  tan"1(lfh>                                   ? 

1-21) 

and  h  is  the  source/detector  distance. 

The  attenuation  coefficient,  u,  used  for  these  calculations  is 

the  "narrow  beam,  total,  without  coherent"  ((— )           )  since 

p  tot,t-coh 

for  calculations  of  crystal  efficiencies  one  is  interested  only  in 
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those  interactions  which  result  in  the  deposition  of  energy.   From 
Fig.  23,  it  is  seen  that  the  intrinsic  efficiency  changes  only 
slightly  for  source  to  detector  distances  greater  than  400  cm.  There- 
fore, Cs-137  and  Co-60  response  functions  were  measured  at  approximately 
400  cm  to  obtain  unfolding  parameters  to  correct  for  this  effect. 
These  measured  response  functions  were  used  to  adjust  the  EF  and  PT 
values  while  leaving  all  other  parameters  unchanged. 

Figure  24  shows  the  unfolded  400-cm  response  functions  while  Table  4 
contains  the  optimized  values  of  the  response  surface  parameters  used 
in  unfolding  the  skyshine  data.  As  with  the  70-cm  unfolded  response 
functions,  the  final  400-cm  results  represent  a  compromise  between 
obtaining  the  known  photopeak  area  and  correctly  compensating  for  the 
Compton  plateau  and  secondary  peaks.   The  analysis  of  the  second  set 
of  response  functions  was  complicated  by  the  presence  of  a  large 
steel  mass  close  to  the  sources  during  the  measurements.   This  resulted 
in  an  abnormally  large  backscatter  peak  which  would  not  appear  in  an 
ideal  response  function.  Therefore,  the  unfolded  400-cm  response  functions 
should  still  contain  a  significant  backscatter  peak  contribution  to 
avoid  over  compensation  when  unfolding  the  skyshine  spectra. 

The  negative  values  observed  in  the  unfolded  results  point  out  a 
weakness  in  the  SEGO  method.   Since  the  stripping  proceeds  from  high 
energy  to  low,  any  errors  in  the  response  matrix  formulation  will  be 
accumulated  as  the  calculation  progresses  and  will  ultimately  be  mani- 
fested in  the  low  energy  portion  of  the  spectrum.  Therefore,  an  accurate 
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Fig.  23.   Intrinsic  efficiency  results  for  the  Nal  crystal  used 
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Table  4.   Optimized  400-cm  Unfolding  Parameters 


E (MeV) 

EF 

PT 

ES 

EC 

BK 

VH 

vu 

VL 

0.0 

1.1810 

1.00 

0.0 

0.0 

0.0 

0.1400 

0.0 

0.0 

0.2 

1.1810 

0.8034 

0.0 

0.0 

0.0120 

0.1000 

0.1738 

0.1050 

0.4 

1.1810 

0.6649 

0.0 

0.0 

0.0120 

0.0600 

0.3480 

0.1200 

0.6 

1.1756 

0.5719 

0.0 

0.0 

0.0130 

0.0240 

0.5481 

0.2000 

0.8 

1.1548 

0.5118 

0.0 

0.0 

0.0140 

0.0225 

0.7352 

0.2000 

1.0 

1.1235 

0.4720 

0.0 

0.0 

0.0143 

0.0234 

0.9224 

0.2200 

1.2 

1.0810 

0.4396 

0.0050 

0.0040 

0.0145 

0.0234 

1.1200 

0.2300 

1.4 

1.0266 

0.4019 

0.0100 

0.0080 

0.0150 

0.0230 

1.2967 

0.2400 

1.6 

0.9597 

0.3464 

0.0150 

0.0120 

0.0150 

0.0230 

1.4838 

0.2500 

Table  5.   Unfolding  Parameters  Used  by  Radiation  Research  Associates" 


E(MeV) 

EF 

PT 

ES 

EC 

BK 

VII 

VU 

VL 

0.0 

1.1750 

1.000 

0.0 

0.0 

0.0 

0.1700 

0.0 

0.000 

0.2 

1.1750 

0.795 

0.0 

0.0 

0.012 

0.1100 

0.164 

0.080 

0.4 

1.1720 

0.680 

0.0 

0.0 

0.012 

0.0750 

0.348 

0.120 

0.6 

1.1550 

0.598 

0.0 

0.0 

0.013 

.0425 

0.540 

0.200 

0.8 

1.1290 

0.530 

0.0 

0.0 

0.013 

0.0355 

0.736 

0.200 

1.0 

1.1020 

0.478 

0.0 

0.0 

0.0135 

0.0225 

0.920 

0.210 

1.2 

1.0800 

0.435 

0.005 

0.004 

0.0145 

0.0229 

1.120 

0.230 

1.4 

1.0610 

0.395 

0.010 

0.008 

0.0145 

0.0229 

1.300 

0.240 

1.6 

1.0440 

0.362 

0.015 

0.012 

0.0150 

0.0229 

1.495 

0.250 
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response  surface  is  imperative  if  an  accurate,  physically  meaningful 
result  is  to  be  obtained,  especially  at  low  energies. 

Since  the  originally  calculated  unfolding  parameters  yielded 
unsatisfactory  results,  a  trial  and  error  approach  was  used  to  adjust  the 
response  surface  parameters  to  obtain  better  results.   Complicating  this 
procedure  is  the  limitation  of  this  method's  accuracy  due  to  its  reliance 
upon  discrete,  idealized  response  functions.   Therefore  the  experimenter 
is  called  upon  to  make  a  judgemental  decision  as  to  whether  or  not  the 
best  possible  result  has  been  obtained  within  the  limitations  of  the 
technique. 

It  has  been  this  author's  experience  that  in  using  the  SEGO  code 
one  must  be  satisfied  with  a  somewhat  less  than  ideal  result.   The 
three  criteria  used  in  adjusting  the  unfolding  parameters  were  proper 
photopeak  area  and  a  peak-to-total  ratio  close  to  unity  without  any 
negative  values  in  the  unfolded  spectra.  While  the  former  criterion 
was  relatively  easy  to  achieve,  it  was  not  possible  to  achieve  the 
latter  criteria.  Thus  an  effort  was  made  to  minimize  the  negative 
contributions  in  the  unfolded  result  although  it  could  not  be  eliminated 
completely. 

Since  there  is  a  great  deal  of  judgement  involved  in  generating 
the  response  matrix  to  be  used  with  the  SEGO  code,  it  is  to  be  expected 
that  two  different  experimenters  working  with  data  from  the  same  spectro- 
meter will  generate  different  sets  of  unfolding  parameters.   As  an 

example,  the  response  matrix  generated  in  this  work  is  not  identical 

3 
to  the  one  published  by  Radiation  Research  Associates  and  which  is  given 
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in  Table  5.  The  magnitude  of  this  effect  is  shown  in  Fig.  25,  where  an 
unfolded  Cs-137  spectrum  is  shown  using  either  the  KSU  or  KRA  parameters. 
While  some  differences  are  apparent,  they  are  certainly  not  significant. 

One  difficulty  encountered  during  the  unfolding  of  the  skyshine 
data  was  the  presence  of  an  iodine  escape  peak  at  approximately  20  keV. 
Since  this  feature  is  a  result  of  a  non-ideal  detector,  it  should  have, 
in  principle,  been  removed  during  the  unfolding.   It  would  have  been 
relatively  easy  to  include  this  structure  in  the  Hyodo  response  functions 
except  that  none  of  the  experimentally  measured  response  functions  con- 
tained an  iodine-escape  peak.   Since  data  was  not  available  from  which 
to  incorporate  the  iodine-escape  peak  into  the  response  matrix,  this 
peak  was  artifically  suppressed  by  linearly  extrapolating  the  spectrum 
to  zero  using  the  two  channels  immediately  above  the  peak  to  define  the 
expected  rapid  decrease  in  the  energy  flux  as  the  spectrum  energy 
decreases  to  zero. 

3.1.5  Determination  of  Gain  and  Zero  Channel 
In  order  to  unfold  the  skyshine  spectra,  the  gain  and  zero- 
intercept  of  each  spectrum  must  be  accurately  known  so  that  the  channel/ 
energy  relationship  may  be  calculated.  A  calibration  spectrum  containing 
five  major  photopeaks  was  taken  before  and  after  each  skyshine  spectral 
measurement.  The  centroid  of  each  photopeak  was  calculated  by  performing 
a  least  squares  fit  of  the  photopeak  to  a  Gaussian.   An  average  spectrum 
gain  and  zero  intercept  was  then  obtained  by  performing  a  linear  least- 
squares  fit  to  the  combined  data  from  both  calibration  spectra.   In  this 
way,  any  gain  shift  that  might  have  occurred  during  the  measurement  due 
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to  an  instrument  or  ambient  temperature  change,  will  be  partially 
compensated  for  by  averaging  the  calibration  data  in  this  manner. 

3.2   Spectral  Comparisons 

One  of  the  primary  purposes  of  the  skyshine  experiment  was  to 
obtain  skyshine  energy  spectra  against  which  calculated  results  could 
be  compared.  With  such  comparisons,  the  accuracy  of  various  photon 
transport  models  can  be  assessed. 

From  August  through  November,  1977,  three  sets  of  spectral  measure- 
ments were  made  of  the  skyshine  radiation  field  produced  by  the  three 
source  configurations  described  in  Chapter  2.   Source/detector  distances 
for  these  measurements  ranged  from  100  to  700  meters  (at  100  m  intervals) 
for  the  unshielded  and  8"  shielded  source  measurements  but  extended  only 
to  500  meters  for  the  16"  shield  cases  before  the  background  signal 
overshadowed  the  skyshine  signal. 

Two  dimensional,  discrete  ordinates  transport  calculations  were 
performed  by  Radiation  Research  Associates  ,  using  the  DOT-3.5  code5 
for  all  source/detector  distances  and  source  configurations.   In  this 
section  the  spectral  results  of  these  calculations  are  compared  to  the 
unfolded  spectra  obtained  with  the  Nal  spectrometer. 

3.2.1  Comparison  of  Unshielded  Spectra 

The  unshielded  source  configuration  consisted  of  an  upward  colli- 
mated  150°  conical  beam  of  Co-60  radiation.   The  unfolded  skyshine 
spectra  (corrected  for  background)  for  this  configuration  are  shown  in 
Fig.  26  for  all  source/detector  distances.   One  immediately  notices  that 
the  spectral  shape  remains  relatively  constant  as  the  source/detector 
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86 
distance  varies  although  some  changes  occurs  as  the  source/detector 
distance  increases.  The  slightly  larger  low  energy  component  in  the 
100  m  spectrum  as  compared  to  the  700  m  result  is  primarily  due  to 
in-silo  scattering  of  the  source  photons.   The  result  of  this  scattering 
would  be  that  the  upward  collimated  radiation  field  would  contain 
multiplely  scattered  low  energy  photons,  as  well  as  uncollided  Ca-60 
gamma  photons.  Each  spectrum  is  also  observed  to  contain  a  rather 
dominant  plateau  region  between  0.25  and  0.75  MeV.  While  being  some- 
what distorted  in  the  100  m  case,  it  is  readily  apparent  in  the  700  m 
example. 

The  effect  of  the  in-silo  scattering  phenomena  on  the  total 
skyshine  exposure  rate  has  been  investigated  by  Ryufuku21  of  the  Japan 
Atomic  Energy  Research  Institute  (JAERI)  using  an  eleven  group  energy 
structure  with  the  DOT-3  transport  code.   The  effect  of  in-silo  scattering 
on  the  skyshine  energy  spectrum  is  shown  in  Fig.  27  for  the  Nal  detector 
geometry.  As  was  expected,  the  inclusion  of  in-silo  scattering  results 
in  a  slight  softening  of  the  spectrum,  particularly  at  small  source/ 
detector  distances.  However,  as  shown  in  Fig.  28,  this  source  of  scattered 
photons  has  a  negligible  effect  upon  the  total  exposure  rate. 

Two  sets  of  DOT  calculations  were  performed  by  RRA  for  the  unshielded 
source  configuration  using  a  two  dimensional  r-z  geometry.   First,  a  39 
energy  group  structure  was  used  with  an  8^  angular  mesh  and  a  P  cross 
section  representation.   Vacuum  boundary  conditions  were  used  except  at 
r=0  where  a  reflection  condition  was  imposed.   To  minimize  streaming  effects 
in  the  calculations,  a  first  collision  source  distribution  was  determined 
and  used  in  the  subsequent  DOT  calculations. 
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Second,  since  the  39  group  calculations  were  extremely  expensive, 
it  was  decided  to  use  a  less  expensive  (and  presumably  less  accurate) 
10  energy  group  model  for  routine  analysis.   In  an  effort  to  preserve 
some  of  the  accuracy  of  the  fine  group  calculations,  a  sensitivity 
study  was  done  in  which  the  original  39  group  energy  structure  was 
collapsed  into  10  groups,  each  of  which  contributed  approximately  the 
same  amount  to  the  total  exposure  rate.   This  optimized  10  group  structure 
was  used  with  a  reduced  angular  grid  and  cross  section  representation 
(S   and  P  ,  respectively)  in  DOT-3.5,  again  with  a  first  collision 
source  distribution.   Complete  details  of  the  modeled  geometry  and 
computational  steps  are  given  in  Reference  3. 

To  aid  in  the  comparison  of  the  39  group  results  with  the  measured 
spectra,  a  Gaussian  smoothing  technique  was  applied  to  the  calculated 
group  fluxes.   With  this  procedure,  the  calculated  spectra  (in  the  form 
of  a  histogram)  can  be  smoothed,  or  smeared,  with  the  same  resolution  as 
the  Nal  spectrometer  which  was  determined  experimentally.   Specifically, 

the  resolution  of  the  Nal  spectrometer  system,  K(E) ,  was  found  to  be 

3 
represented  by 

K(E)  =  0.0745E"0-44  (3.2-1) 

where  R(E)  is  the  fractional  resolution,  or  line  width,  defined  as 

K(E)  =  BK  (3.2-2) 

max 

where  C    is  the  pulse  height  (MeV)  corresponding  to  the  centroid  of 

the  photopeak  and  FWHM  is  the  pulse  height  width  (MeV)  between  the 
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points  at  which  one-half  the  maximum  value  occurs  (i.e.,  the  full 
width  at  half  maximum) .   Since  the  photopeak  for  a  Nal  spectrum  can 

be  well  represented  by  a  Gaussian,  the  quantity  K  can  be  related  to 

19 
the  standard  deviation  of  a  normal  distribution  with  mean  E  by 


2.35  o(E) 
KOJj) g- *-  (3.2-3) 

Therefore,  denoting  the  smoothed  differential  flux  in  channel  i  by  $  , 
and  the  calculated  (histogram)  differential  flux  in  group  j  by  $.,  the 
smoothed  spectrum  can  be  represented  by 
NG 


♦  =  I  S   *       1-1,2... .8 

j-1  3  J 


(3.2-4) 


where  NG  is  the  number  of  group  fluxes  in  the  calculated  spectra,  N  is 

the  number  of  channels  in  the  smeared  spectrum,  and 

AE. 
1 

1J   a,  /2T  "    "  *  *' 


S,,  -  *-  exp  [-(E,-E,)2/2a.2]     .        (3.2-5) 


The  result  of  this  smearing  procedure  on  the  calculated  group  fluxes 
is  a  "smooth  spectrum"  which  should  very  closely  approximate  the  unfolded 
results  if  the  transport  model  is  sufficiently  accurate.   Unfortunately 
this  procedure  could  not  be  applied  to  the  10  group  results  because  the 
energy  group  width  used  in  the  calculations  was  much  larger  than  the 
resolution  width  of  the  spectrometer  and  such  a  smoothing  technique  would 
introduce  artificial  structure  into  the  smoothed  spectrum. 


91 

A  comparison  between  the  calculated  and  unfolded  incident  un- 
shielded skyshine  energy  spectra  is  shown  in  Fig.  29  for  three  repre- 
sentative source/detector  distances.   The  smooth  curves  with  identifying 
symbols  are  the  unfolded  Nal  results.   (The  symbols  are  located  at  every 
fifteenth  data  point.)  The  other  smooth  curves  (without  symbols)  are  the 
39  group  Gaussian  smoothed  results  calculated  by  DOT,  while  the  histograms 
represent  the  10  group  results.   In  the  100  meter  case,  the  agreement  is 
very  good  for  energies  above  0.50  MeV.   Below  this  point  the  measured 
spectrum  is  considerably  softer  than  that  of  the  DOT  results.  This 
difference  was  expected  since  in-silo  scattering  was  not  included  in  the 
DOT  calculations  done  by  ERA.   The  10  group  calculation  agrees  fairly 
well  but  slightly  overestimates  the  measured  result  in  the  first  few 
groups  and  is  also  deficient  in  describing  the  70  keV  characteristic 
skyshine  peak.   Both  of  these  effects  are  apparently  attributable  to 
the  energy  bins  being  too  wide  to  closely  approximate  the  actual 
incident  spectrum. 

The  agreement  between  the  predicted  and  measured  spectra  is 
excellent  in  the  400  meter  case.  The  10  group  result  again  slightly 
overestimates  at  high  energies  and  also  lacks  sufficient  detail  in  pre- 
dicting the  low  energy  peak.  The  Gaussian  smoothed,  39  group,  calculated 
result  is  almost  indistinguishable  from  the  measured  spectra.   The  in-silo 
scattering  effect  should  be  negligible  at  this  distance  since  the  lower 
energy,  silo-scattered  photons  are  preferentially  removed  from  the  sky- 
shine radiation  at  large  source/detector  distances. 
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Fig.  29.   Comparison  of  representative  calculated  and  unfolded 

measured  spectra  for  the  unshielded  source  configuration. 
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The  agreement  in  the  700  meter  results  is  somewhat  difficult  to 
interpret  since  the  source/detector  distance  at  which  the  measurement 
was  made  is  approximately  5%  less  than  the  distance  used  in  the 
corresponding  DOT  calculation.  Disregarding  the  differences  in  distances, 
the  results  appear  to  agree  very  well. 

3.2.2  Comparison  of  Shielded  Spectra 

The  shielded  source  configurations  consisted  of  a  Co-60  source 
covered  first  by  8"  and  then  16"  of  concrete.   Only  ten  group,  DOT 
calculations  were  performed  for  the  shielded  cases.  Again  an  S., 
angular  grid  and  P  cross  section  representation  were  used.  To  reduce 
further  the  computation  time,  the  calculations  were  broken  into  several 
stages  by  RRA.   First,  the  radiation  penetrating  the  source  shield  was 

calculated  with  the  one  dimensional,  discrete  ordinates  transport  code, 

22 
ANISN.     This  shield  penetration  radiation  was  then  used  as  a  source  for 

an  air  skyshine  calculation  using  DOT. 

The  unfolded  results  for  the  8"  and  16"  shielded  source  config- 
urations are  shown  in  Figs.  30  and  31,  respectively.   As  with  the 
unshielded  results,  each  spectrum  peaks  at  approximately  75  keV,  but 
the  middle  energy  range  plateau  is  no  longer  present.   Instead  the 
unfolded  spectra  drop  off  in  an  almost  perfect  exponential  fashion 
above  the  low  energy  peak  resulting  in  somewhat  softer  spectra  than 
in  the  unshielded  case.   Although  the  100  meter  spectra  are  still 
slightly  distorted  by  the  in-silo  scattered  component,  this  scattered 
source  effect  is  considerably  decreased  by  the  silo  roof  shields 
compared  to  the  unshielded  skyshine  results. 
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Fig.  30,   Comparison  of  unfolded  Nal  spectra  for  the  8" 
hielded  source  configuration. 
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Fig.  31.  Comparison  of  unfolded  Nal  spectra  for  the  16"  shielded 
source  configuration. 


96 


In  Fig.  32  an  unfolded  skyshine  spectrum  measured  at  100  meters 
with  the  16"  shielded  source  configuration  is  shown.   The  feature  to  note 
in  this  spectrum  is  the  plateau  occurring  between  0.80  and  1.30  MeV, 

presumably  caused  by  direct  penetration  of  the  silo  wall  by  source  photons. 

21 
This  effect  was  also  studied  by  Ryufuku   to  determine  its  contribution 

to  the  measured  exposure  rate.   Again  using  the  DOT-3  transport  code,  it 

was  found  that  for  source/detector  distances  less  than  five  meters,  the 

direct  penetration  component  accounts  for  approximately  15%  of  the  total 

exposure  rate  while  at  eighteen  meters  this  contribution  falls  below  1%. 

These  results  for  the  wall-penetration  component  are  shown  in  Fig.  33 

along  with  some  experimental  results  obtained  with  a  high  pressure 

ionization  chamber  (which  is  described  in  detail  in  Chapter  4) .  The  large 

differences  between  these  calculated  and  measured  values  probably  arises 

because  of  the  sensitivity  of  the  calculated  results  to  the  assumed 

problem  geometry  which  was  an  approximation  of  the  actual  source-silo 

configuration.   The  effect  of  this  geometrical  approximation  would 

undoubtedly  be  most  pronounced  at  small  source/detector  distances. 

Nevertheless,  these  results  indicate  that  the  silo  wall  penetration 

component  is  undoubtedly  small  compared  to  the  skyshine  component. 

Representative  spectra  for  the  8"  shielded  source  configuration  are 

shown  in  Fig.  34.  Again  it  is  seen  that  the  10  group  calculations 

seriously  underpredict  the  70  keV  peak  in  the  measured  spectra.   The 

apparent  close  agreement  for  intermediate  and  high  energies  is  somewhat 

misleading  because  of  an  error  in  the  DOT  modeling  of  the  roof  shield  in 

the  8"  shield  configuration.   The  actual  thickness  of  the  shield  in  this 
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g.  32.  Unfolded  Nal  spectrum  showing  direct  penetration  plateau 
for  the  16"  shielded  source  configuration.   Source  to  de- 
tector distance  is  100  meters. 
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Fig.  33.   Direct  penetration  exposure  rate  contribution  for  various 
detector  heights,  h,  above  grade. 
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Fig.  34.   Comparison  of  representative  calculated  and  unfolded 
measured  spectra  for  the  8"  shielded  source  con- 
figuration. 
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configuration  was  21.03  cm  whereas  20.03  cm  was  used  in  the  DOT  calcula- 
tions.  By  performing  a  simple  two-point  exponential  fit  to  the  exposure 
rates,  it  was  estimated  that  this  error  causes  an  overprediction  of  10  to 
20%.   Other  estimates  place  the  error  at  15%  . 

Finally,  Fig.  35  shows  a  comparison  of  representative  16"  shielded 
results.  As  expected,  the  agreement  is  very  good  for  all  energy  regions 
except  below  100  keV. 

It  should  be  noted  that  for  all  source  configurations,  the  shape 
of  the  DOT  calculated  spectra  are  extremely  good,  especially  in  the 
39  group  results.  However,  the  normalization,  or  magnitude,  of  these 
spectra  is  inadequate  for  the  10  group  calculated  spectra.  Additionally, 
negative  flux  values  were  obtained  in  the  first  few  energy  groups  of  the 
calculated  results,  which  were  suppressed  to  zero  in  this  analysis. 
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Comparison  of  representative  calculated  and  unfolded 
measured  spectra  for  the  16"  shielded  source  con- 
figuration. 
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4.0  Measurement  and  Comparison  of  Exposure  Rates 

In  addition  to  the  skyshine  spectral  measurements  which  served 
as  a  test  of  the  various  photon  transport  models,  the  second  major  ob- 
jective of  the  skyshine  experimental  program  was  the  measurement  of  the 
total  skyshine  exposure  rate  as  a  function  of  source/detector  distance 
for  both  the  shielded  and  unshielded  3ource  configurations.   Such  exposure 
rate  measurements  serve  as  additional  benchmark  data  against  which  various 
calculational  schemes  for  the  prediction  of  the  skyshine  exposure  rates  can 
be  compared.   These  measurements  may  also  be  used  directly  as  design  data 
for  facilities  in  which  skyshine  radiation  is  of  concern. 

In  this  phase  of  the  experiment,  a  high  pressure  ionization 
chamber  (HPIC)  was  used  to  measure  the  Atr  skyshine  exposure  rate  (I.e., 
omnidirectional  exposure  rate)  over  a  700  meter  range  of  source-to- 
detector  distances.   The  source  configurations  employed  are  the  same 

93 
as  those  used  for  the  skyshine  spectral  measurements""  .  A  10  inch 

(25.4  cm)  diameter,  spherical,  twenty-five  atmosphere  (2.5  MPs)  argon- 
filled  ionization  chamber  (Reuter-Stokes  Area  Monitor  System,  Model 

24  25 
RS-111)      was  used  to  measure  the  exposure  rate  resulting  from  the 

skyshine  radiation  field.   Results  and  details  of  these  measurements 

3  23 
are  given  elsewhere  ' 

The  exposure  rate  measured  with  the  Reuter-Stokes  detector  must 

be  corrected  to  obtain  the  true  exposure  rate  because  the  system 

response  is  not  linear  with  energy,  especially  below  300  key.  Although 
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the  system  response  curve  was  supplied  by  the  manufacturer,  the 
incident  energy  spectrum  is  also  required  to  calculate  a  correction 
factor  which  compensates  for  the  non-ideal  behavior  of  the  system. 
Initial  tests  indicated  that  the  ionization  chamber  was  insensitive 
to  the  direction  of  the  incident  radiation  and,  therefore,  only  the 
energy  dependence  of  the  scalar  flux  is  necessary  for  the  calculation 
of  this  correction  factor. 

While  the  Nal  spectral  measurements  provided  some  information 
about  the  energy  spectrum  of  the  incident  radiation  field,  the  Nal 
detector  was  collimated  and  hence  the  energy  spectrum  of  only  a  portion 
of  the  total  radiation  field  was  obtained.   Consequently,  the  original 

correction  of  the  exposure  rate  data  by  RRA  for  non-ideal  response  of  the 

3 

HPIC  was  based  on  the  DOT-3.5  calculated  4n  energy  spectra.   In  this 

chapter,  the  calculation  of  the  spectrum  correction  factors  for  the 
HPIC  results  is  re-examined  and  compared  to  a  correction  procedure 
which  uses  Nal  spectral  measurements,   the  Nal  spectral  measurements 
are  also  used  to  predict  an  exposure  rate  for  the  4ir  skyshine  radiation 
field.   These  collapsed  Nal  measurements  then  serve  as  an  independent 
determination  of  the  total  exposure  rates  against  which  the  corrected 
Reuter-Stokes  HPIC  values  may  be  compared. 

4.1  Correction  of  Observed  Exposure  Rates 
4.1.1  Theoretical  Considerations 
In  general,  the  exposure,  X,  for  a  given  photon  fluence  is 
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defined  as 


M 


(4.1-1) 


where  AQ  is  the  sum  of  the  electrical  charges  on  all  the  ions  of  one 
sign  produced  in  air  when  all  the  electrons  (negatrons  and  positrons) 
liberated  by  photons  in  a  volume  element  of  air,  whose  mass  is  £m, 
are  completely  stopped  in  air.  The  unit  of  exposure  is  the  roentgen 
(R)  which  is  equivalent  to  2.58  x  10~  Coulombs/kg.  It  follows  that 
the  exposure  rate  is  simply  the  amount  of  exposure  per  unit  time  and 
may  be  calculated  as 

•    _  E*e   (u   /P)alr  ,         (4.1-2) 

X  =  — "r    en 

where  E  is  the  energy  of  the  incident  photons,  4>  is  the  photon  flux, 
e  is  the  electron  charge  (1.6  x  10    C) ,  W    is  the  average  amount 

of  energy  required  to  create  one  ion-electron  pair  (equal  to  approx- 

air 
imately  34  eV  over  the  photon  energy  range  of  interest),  and  (u  /p) 

is  the  mass  energy  absorption  coefficient  for  air.   The  use  of  this 

particular  coefficient  implies  the  loss  of  all  secondary  photons,  e.g., 

bremsstrahlung.   The  substitution  of  numerical  values  into  Eq.  (4.1-2)  yields 

X  =  65.664  E  *(u  /p)air  ,    (4.1-3) 

en 

where  X  is  in  yR/h,  E  in  MeV,  $  in  cm   sec   and  (y  /p)    in  cm  g  . 
When  the  incident  spectrum  is  represented  by  discrete  components  (e.g., 
a  pulse  height  spectrum),  Eq.  (4.1-3)  can  be  formulated  as 
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NC 
X  =  65.664  I     E  ♦.  (u  /p)?lr    ,    (4.1-4) 

1=1     x   en   x 

where  NC  represents  Che  number  of  components  in  the  spectrum,  E.  and 

air 
(wen/p)i   are  appropriately  averaged  values  for  the  i-th  energy  interval 

in  the  spectrum,  and  $  is  the  energy-integrated  flux  over  the  i-th  interval. 

In  an  ionization  chamber  such  as  the  one  used  in  the  skyshine 
experiment,  the  exposure  rate  is  Inferred  from  the  current  produced 
by  the  radiation-ionized  filling  gas.   In  the  ideal  case,  the  response 
of  the  detector  would  be  directly  related  to  the  radiation  intensity 
by  some  conversion  constant,  regardless  of  the  energy  of  the  ionizing 
radiation.   In  practice,  the  response  of  the  system  is  not  ideal  because 
the  steel  sphere  filled  with  pressurized  argon  does  not  respond  as 
would  a  detector  composed  totally  of  air,  i.e.,  it  is  not  "air  equiva- 
lent".  Therefore,  in  general,  an  energy  dependence  correction  factor 
must  be  derived  to  correct  the  observed  exposure  rate  readings. 

The  current  collected  from  the  ionized  filling  gas  1(E),  in  amperes, 
due  to  an  incident  radiation  field,  $(E),  is  given  by 

1(E)  =  K(E)-X(E)     ,  (4.1-5) 

where  X(E)  is  the  exposure  rate  as  a  function  of  photon  energy  in 
pR/h  and  K(E)  is  a  conversion  constant  which  can  be  expressed  as 

K(E)  -  CF(E)  .  (4.1-6) 

In  the  above  equation,  C  is  a  proportionality  constant  determined 
at  some  specific  reference  energy  and  F(E)  is  the  response  of  the  argon- 
filled  sphere  with  energy  compared  to  the  reference.   The  measured 
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energy  sensitivity  of  the  system,  F(E),  is  shown  in  Fig.  36,  for  the 

chamber  size  and  filling  pressure  used  in  this  experiment. 

From  Eq.  (4.1-5),  the  total  current  collected,  I  ,  due  to  an  in- 

Y 
cident  polyenergetic  radiation  field  is 

E  E 

rmax  max 

I  -     dE  1(E)  =     dE  K(E)X(E)  ,     (4.1-7) 

'   ;0  ->0 

Where  Emax  is  the  maximum  photon  energy  in  the  spectrum.   After 
substitution  of  Eq.  (4.1-6),  I  can  be  written  as 

E 
(•max 

Iy  -  C  I    dE  F(E)X(E)  .  (4.1-7) 

To  obtain  an  expression  for  X(E),  Eq.  (4.1-3)  may  be  expressed  in 
terms  of  continuous  functions  as 


X(E)  =  65.664  {E<j>(E)(u   (E)/p)}  ,       (4.1-8) 


and  upon  substitution  in  Eq.  (4.1-7),  I  becomes 

Y 
E 
max 


I  -  65.664C      dE  F(E)E*(E)(u   (E)/p).        (4.1-9) 
'  '0  en 

For  convenience,  an  average  system  sensitivity,  F,  may  be  defined  as 

E 
r   max 

dE  F(E)E<|>(E)<u   (E)/p) 

J  =  lo en 

F  "  ~i (4.1-10) 

r   max 

dE  E<KE)(y   (E)/p) 
In  en 
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which  can  be  physically  interpreted  as  the  ratio  of  the  observed 
exposure  rate  with  a  non-ideal  system  sensitivity  to  the  exposure 
rate  obtained  with  an  ideal  system  sensitivity.  With  this  definition, 
Eq.  (4.1-9)  becomes 

\   =  CFX       ,  (4.1-11) 


where 


E 
/•max 

Jo 


dE  E  ij>(E)  (u   (E)/p)   .         (4.1-12) 


Since  the  argon-filled  sphere  is  surrounded  by  a  cubical  aluminum 
protective  housing  which  was  not  present  in  the  original  determination 
of  F(E),  an  additional  minor  correction  needs  to  be  made  to  account 
for  the  attenuation  of  the  incident  radiation  by  this  housing.   There- 
fore, Eq.  (4.1-10)  is  modified  to 

f  air  "'Vi*  C 

F  =  ii^ ' ,  (4.1-13) 

I   E,  *,  (y   /p)fir 
.~.    i  i   en   x 

where  F  has  been  expressed  in  terms  of  a  discrete  energy  representation 
and  t  is  the  wall  thickness  of  the  aluminum  shell. 

An  expression  for  the  exposure  rate,  X,  can  now  be  obtained  from 
Eq.  (4.1-11)  as 

I. 

CF 
For  a  specific  reference  or  calibration  energy,  denoted  by  c,  Eq. 
(4.1-6)  becomes 


X  -  -*  •  (4.1-14) 
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^c  "' 

K  -  CF  e  (4.1-15) 


where  the  effect  of  the  attenuation  by  the  aluminum  housing  has 

been  included.   Thus 

-(y  )   -t 
K  e 


(4.1-16) 


Substitution  of  Eq.  (4.1-16)  into  Eq.  (4.1-14)  then  yields 

.   I  F  e  a  c 

X  =  -JL£ .  (4.1-17) 

K  F 

c 

Furthermore,  it  should  be  noted  that  the  current  measured  by 

the  ion  chamber  is  the  sum  of  two  components.   The  first,  I  ,  is 

due  to  the  ionization  of  the  argon  filling  gas  by  the  incident  radiation 

field.   The  second,  I  ,  is  a  result  of  inherent  svstem  noise  and  back- 
b 

ground  radiation.  Thus 


IY  "  It  "  \  (4-1-18) 

where  I  is  the  total  current  collected.   Substitution  of  this  ex- 
pression into  Eq.  (4.1-17)  yields 

.         (I  -  I,  )  F  e  a  C 

X 2— 3s .  (4.1-19) 

K  F 
c 

From  this  result,  it  is  observed  that  the  true  exposure  rate 

as  given  by  Eq.  (4.1-19)  is  composed  of  two  parts,  the  first  being 

the  observed  or  measured  exposure  rate, 
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c 
and  the  second  being  a  correction  factor  to  account  for  the  non-ideal 
response  of  the  system  with  energy.  Equivalently, 

•   -  <rr  "  V 

X  ■  1     R  (4.1-21) 

c 

where  f  is  a  "system  energy  correction  factor"  defined  as 

F  e 
?-— .  (4.1-22) 


This  correction  factor,  f,  is  thus  used  to  correct  the  HPIC  reading  to 

yield  the  true  exposure  rate.  Notice  that  to  obtain  a  value  for 

f  (or  F)  ,  both  the  system  energy  sensitivity,  F(E),  and  the  incident 

energy  spectra,  t|>(E),  must  be  known.   The  effect  of  this  non- ideal 

behavior  on  the  exposure  rate  is  shown  in  Figs.  37-39  for  the  unshielded, 

8"  shielded  and  16"  shielded  skyshine  source  configurations,  respectively. 

The  quantities  shown  on  these  graphs  are  the  exposure  rate  contributions 

for  measured  Nal  energy-dependent  fluxes  (unfolded)  as  computed  for 

both  the  ideal  (i.e.,  F(E)=1)  and  the  non-ideal  system  energy  sensitivities. 

The  non-ideal,  or  Reuter-Stokes  HPIC  exposure  rate  contribution  from  the 

flux  in  energy  interval  i  is  defined  as 

Xs1  ,_a    "   65.66    I*,.,    E.j    4.^,     (u      /p),,,  (4.1-23) 

contr.i+Sj  {   i+Jj     tf]  Yi+Ss   v  en       i+ijj  v  ' 

while   the  ideal  exposure   rate   contribution   is   calculated  with  F. ,,  =1 
for  all  i.      The   subscript   "i+V  denotes    the  value  of   the   function 
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evaluated  at  the  center  of  the  energy  interval  which  for  these  figures 
had  a  width  of  5  keV.   In  these  figures,  the  area  under  the  ideal 
contribution  curve  is  equal  to  the  true  exposure  rate  as  would  be  calcu- 
lated from  the  unfolded  Nal  spectrum,  while  the  area  under  the  non-ideal 
response  curve  is  equal  to  the  exposure  rate  which  would  have  been 
measured  by  the  Reuter-Stokes  HPIC  detector,  if  it  had  experienced  the 
same  degree  of  collimation  as  did  the  Nal  detector. 

Because  of  the  peaked  nature  of  F(E)  at  approximately  100  keV,  the 
Reuter-Stokes  detector  overemphasizes  the  low  energy  contribution  to 
the  exposure  rate.  As  opposed  to  this,  the  middle  and  high  energy  con- 
tributions are  underestimated  because  the  value  of  F(E)  is  less  than 
unity  in  these  regions  (.88  at  1.25  MeV) .  The  extreme  low  energy  con- 
tributions (E  <  50  keV)  is  completely  lost  to  the  Reuter-Stokes  chamber 
because  of  the  cutoff  in  the  sensitivity  at  this  energy.   Since,  as  dis- 
cussed in  Chapter  3,  the  skyshine  spectra  peaked  at  approximately  70  keV, 
the  ionization  chamber  may  miss  a  significant  portion  of  the  incident 
spectrum.  However,  this  cutoff  in  the  detector  sensitivity  did  not 
turn  out  to  be  serious  since,  as  can  be  seen  in  Figs.  37-39,  the  maximum 
contribution  to  the  exposure  rate  comes  at  600  and  300  keV  for  the  un- 
shielded and  shielded  skyshine  source  configurations,  respectively. 

4.1.2  Calculation  of  the  Correction  Factor  f 

Two  problems  are  encountered  in  the  calculation  of  the  HPIC  correction 
factor  f.   The  first  problem  is  that  the  true  incident  spectrum,  i)i(E), 
or  <ji   is  generally  not  known.   Secondly,  large  errors  may  be  incurred 
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by  Che  representation  of  the  rapidly  varying  continuous  quantities 
in  Eq.  (4.1-7)  by  the  discrete  or  energy-group  averaged  values  in 
Eq.  (4.1-13). 

The  simplest  method  of  obtaining  group  values  for  the  discrete 
representation  is  to  evaluate  the  continuous  function  at  the  midpoint 
of  the  energy  bin  and  thereafter  assume  that  it  is  constant  for  that 
particular  energy  range.   Such  a  simple  discretization  procedure  may 
be  acceptable  if  the  energy  mesh  or  bin  width  is  sufficiently  fine. 
A  more  sophisticated  and  accurate  way  of  treating  this  problem  is  to 
use  spectral  weighting  to  average  a  particular  parameter  across  each 
energy  bin.  For  example,  given  some  energy  dependent  quantity  B(E), 
the  group  averaged  value  is  given  by 

r%  +ae 

g    g 
J    dE  8(E)$(E) 
E 
a     =  _S 


E  +AE 
g   g 

dE  *(E) 


(4.1-24) 


g 

where  E  in  the  present  application  corresponds  to  the  lower  edge  of 
some  energy  bin  "g"  and  AE  is  the  width  of  that  bin. 

Just  as  in  Eq.  (4.1-10),  Eq.  (4.1-24)  involves  the  incident 
spectrum,  § (E) ,  which  is  generally  unknown  and  therefore  an  approximate 
representation  must  often  be  used.   The  philosophy  of  spectral  weighting 
involves  the  use  of  a  "hyper-fine"  discrete  representation  of  if(E)  to 
calculate  group  constants  to  be  used  with  coarse  energy  group  represen- 
tations of  the  incident  spectrum.   The  bin  widths  in  the  hyper-fine 
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representation  need  to  be  sufficiently  small  so  that  the  variation  of 
3(E)  can  be  considered  negligible  across  any  one  bin  and  therefore  mid- 
point values  can  be  used. 

Two  hyper-fine  mesh  approximations  to  4>(E)  were  available  in  this 
analysis,  the  first  being  the  Gaussian  smoothed  39  group  DOT  results 
and  the  second  being  the  unfolded  Nal  spectra.  The  Gaussian  smoothed 
calculated  spectra  represented  the  best  estimates  available  of  the  in- 
cident flux  for  the  unshielded  source  configuration.  As  discussed  in 
Chapter  3,  the  39  group  spectra  agreed  extremely  well  with  the  unfolded 
Nal  results  and  thus  one  assumes  that  the  prediction  of  the  4ir  geometry 
incident  flux  would  be  equally  good  although  experimental  data  were  not 
available  for  a  direct  comparison. 

However,  unfolded  Nal  spectra  represented  the  only  hyper-fine  mesh 
estimate  of  the  incident  flux  for  the  shielded  source  configurations, 
since  only  coarse  mesh,  10  group  DOT  calculations  were  performed  for 
these  cases.  Because  the  4it  energy  spectra  were  expected  to  be  "softer" 
than  the  collimated  geometry  spectra  (an  observation  from  the  39  group 
DOT  results),  some  error  may  be  incurred  by  approximating  the  4tt 
incident  flux  by  a  collimated  result. 

The  calculation  of  f  for  the  unshielded  source  configuration  pro- 
vided an  excellent  opportunity  to  determine  the  sensitivity  of  the  group 
constants  to  the  assumed  spectral  shape  used  for  weighting  since  the 
values  obtained  using  the  collimated  Nal  spectra  as  weighting  functions 
could  be  compared  to  those  obtained  using  the  Gaussian  smoothed  39  group 
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results  as  the  weighting  spectra.   The  unshielded  calculations  also 
provided  an  excellent  determination  of  the  sensitivity  of  the  result 
to  the  coarse  group  approximation  of  $(E)   used  in  Eq.  (4.1-13).  With 
the  39  group  results  as  a  benchmark,  the  results  obtained  with  the  10 
group  approximation  could  thus  be  examined.   It  was  imperative  that  the 
sensitivity  of  the  calculated  f's  to  these  parameters  be  determined 
with  the  calculated  unshielded  spectra  because  the  39  group  calculations, 
which  serve  as  the  best  estimate  of  the  4n  spectrum,  were  not  performed 
for  the  shielded  source  configurations. 

The  procedure  followed  in  the  determination  of  the  group  constants 
required  for  the  calculation  of  f  was  to  develop  discrete  hyper-fine 
mesh  representations  of  F,  (u  /p)a  r,  (u  /p)A  ,  and  $.  The  represen- 
tations of  cji  have  already  been  discussed.  It  was  assumed  that  the  bin 
width  of  the  hyper-fine  mesh  structure  was  sufficiently  narrow  so  that 
negligible  error  was  incurred  by  approximating  continuous  functions  by 
their  value  at  the  centroid  of  each  bin,  i.e., 

E  +AE 
g       g 

dE  S(EJ*<B) 

E 

jf^SF =    S(Eg+AEg/2)      ,  (4.1-25) 

J  dE*(E) 

E 
g 

for  sufficiently  small  AE  . 

The  data  used  in  the  construction  of  the  hyper-fine  mesh  repre- 
sentation of  the  photon  coefficients  were  taken  from  Reference  20.  A 
two-point  exponential  fit  was  used  successively  between  these  data 
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points  to  obtain  values  at  5  keV  intervals  corresponding  to  the  mid- 
points of  the  energy  bins  in  the  hyper-fine  structure. 

The  discrete  representation  of  F(E)  was  accomplished  through  the 
use  of  data  from  Reference  25  which  are  tabulated  in  Table  6.  As  des- 
cribed above,  a  two-point  exponential  fit  was  used  successively  between 
the  data  points  in  the  energy  regions  from  50  to  80  keV  and  from 
0.15  to  2.0  MeV.  A  three-point  quadratic  fit  was  used  between  80  and 
150  keV  to  approximate  the  peaked  portion  of  F(E).   The  value  of  Fc  and 
GO   were  taken  at  the  average  Co-60  emission  energy  of  1.25  MeV. 

Table  6.   Data  Used  in  Construction  of  F(E) 


Energy  (keV) 

F(E) 

50 

0.064 

60 

0.408 

80 

1.390 

100 

1.562 

150 

1.270 

300 

0.982 

500 

0.940 

1000 

0.886 

2000 

0.886 

Since  the  data  in  this  analysis  consisted  of  discrete  components, 
Eq.  (4.1-24)  must  be  written  as 
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NB 

I  * 

before  performing  the  spectral  weighting.   In  Eq.  (4.1-26),  NB 
stands  for  the  number  of  hyper-fine  bins  in  energy  group  g.  As  can 
be  seen  from  Eq.  (4.1-13),  two  sets  of  group  constants  are  necessary 
for  the  calculation  of  f.  With  reference  to  Eq.  (4.1-26),  the  first 
set  of  group  constants  corresponds  to 

-OOi  ' 
B.  -  F.  E.  (u  /p)  .  e         ,  (4.1-27) 

i    11   en    i  ■  i^»*  •»* 

while  the  second  set  corresponds  to 

3i  =  Ei  (wen/p)i  "  (4-1"28) 

Selected  results  of  this  sensitivity  study  for  the  unshielded 
skyshine  case  are  presented  in  Table  7.   The  quantities  tabulated 
are  the  f  correction  factor  calculated  from  Eq.  (4.1-18)  as  a  function 
of  source/detector  distance  and  weighting  spectrum.  Results  are  also 
shown  comparing  the  values  obtained  when  the  10  group  approximation 
to  $(E)  is  used  in  Eq.  (4.1-8)  instead  of  the  fine  mesh  39  group 
spectra.  The  second  subscript  on  f  refers  to  the  weighting  spectrum 
used  while  the  superscript  and  the  first  subscript  refer  to  the  approxi- 
mation to  $(B)  used  in  Eq.  (4.1-8).   For  example,  !*?  1QQmGS   signifies  an 
f  calculated  using  a  4n  geometry,  10  group  DOT  calculated  approximation 
to  4>(E)  with  group  constants  computed  from  a  100  meter  Gaussian  smoothed 
weighting  spectrum. 
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As  can  be  seen  from  the  results,  Che  variations  in  the  f  values 
are  negligible  regardless  of  whether  a  4ir  or  collimated  geometry 
weighting  spectrum  is  used.   This  is  somewhat  surprising  since  the 
shapes  of  the  two  spectra  are  significantly  different.   It  is  also 
observed  that  the  value  of  f  changes  only  slightly  when  a  10  group 
spectrum,  with  properly  average  group  constants,  is  used  in  Eq.  (4.1-8) 
instead  of  the  39  group  approximation.   The  f  values  shown  in  Table  8 
are  taken  from  Reference  3  for  the  unshielded  source  configuration. 
The  differences  between  the  values  produced  in  this  work,  and  those  shown 
in  Table  8  are  less  than  5%.   The  calculational  steps  used  in  obtaining 
the  values  in  Table  8  are  not  described  in  Reference  3  and  therefore, 
it  is  not  known  why  these  apparent  differences  arise.   There  is  certainly 
room  for  error  in  the  construction  of  F(E)  due  to  the  sparseness  of  data 
supplied  in  Reference  25. 

The  close  agreement  of  the  results  contained  in  Table  1  suggests 
that  for  the  analysis  of  skyshine  data  one  may  use  collimated  unfolded  Nal 
spectra  in  the  calculation  of  group  constants  when  a  4tr  weighting  spectrum 
is  unavailable.   Since  39  group  calculations  were  done  only  for  the  un- 
shielded source  configuration,  unfolded  Nal  spectra  must  be  used  as 
weighting  spectra  in  the  calculation  of  the  group  constants  for  the 
shielded  source  configurations.   Additionally,  the  relatively  good  agree- 
ment between  the  10  and  39  group  calculated  ?'s  for  the  unshielded  case 
gives  one  a  certain  amount  of  confidence  in  the  results  obtained  with  the 
same  technique  for  the  shielded  source  configurations. 
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Table  8. 

Published  Values  of  Che  Energy  Dependence  Correction 
Factor  for  the  Unshielded  Source  Configuration 

Radial                 j 
Distance  (m) 

50                  .885 

100                 .905 

200                 .901 

300                 .895 

400                 .890 

500                  .887 

600                 .883 
700                 .883 
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4.1.3  Correction  Factor  Values  for  the  Unshielded  Source  Configuration 
To  correct  the  HPIC  measurements  for  the  non-ideal  energy  response 
of  the  system,  f  values  were  calculated  at  100  meter  intervals  along  the 
700  meter  measurement  baseline.  The  weighting  spectra  used  were  the 
Gaussian  smoothed,  39  group,  4ir  geometry  DOT  results.   The  results  of 
these  calculations  are  given  in  Table  9  where  they  are  also  compared  to 
those  given  in  Reference  3.   From  these  restuls  it  is  observed  that  the 
f  values  obtained  in  this  work  increase  slightly  with  increasing  source/ 
detector  distance.   There  are  two  competing  processes  which  contribute 
to  this  trend  which  can  be  seen  in  Fig.  37.  The  peak  portion  (70-200  keV) 
of  the  non-ideal  energy  sensitivity  curve  tends  to  force  f  below  unity 
since  the  exposure  rate  as  seen  by  the  HPIC  is  larger  than  reality.  How- 
ever, the  opposite  effect  occurs  above  approximately  250  keV  where  the 
true  exposure  rate  is  underestimated  by  the  HPIC  thus  tending  to  make  f 
greater  than  unity. 

Both  of  these  processes  are  significant  in  the  unshielded  source  con- 
figuration results.  As  can  be  observed  in  Fig.  26,  the  number  of  photons 
in  the  region  immediately  above  the  75  keV  maximum  decreases  with  increasing 
source/detector  distance  thus  tending  to  make  f  closer  to  unity.  The  ratio 
of  the  number  of  photons  at  75  keV  to  the  number  of  photons  at  750  keV  is 
increasing  with  source/detector  distance  which  tends  to  make  f  less  than 
unity.   The  former  effect  dominates  in  these  bare  skyshine  spectra  which 
explains  the  observed  upward  trend  in  the  f  values  (although  the  latter 
effect  tends  to  dampen  this  variation) .   It  is  interesting  to  note  that 
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the  RRA  calculated  values  exhibit  the  opposite  behavior,  a  trend  which 
seems  to  have  little  physical  justification. 

Table  9.   Comparison  of  Unshielded  f   Values 


Radial  Distance 

This 

3 
RRA 

(meters) 

Work 

100 

.881 

.905 

200 

.892 

.901 

300 

.902 

.895 

400 

.909 

.890 

500 

.913 

.887 

600 

.916 

.883 

700 

.918 

.883 

4.1.4  Correction  Factor  Values  for  the  Shielded  Source  Con- 
figurations 

Unfolded  Nal  spectra  were  used  as  the  weighting  spectra  in  calcu- 
lating the  f's  for  the  shielded  source  configurations.   As  in  the  un- 
shielded case,  these  values  were  calculated  at  100  meter  intervals 
along  the  700  meter  measurement  baseline.  The  results  of  these 
calculations  for  the  8"  shielded  source  configuration  are  shown  in 
Table  10  as  well  as  the  corresponding  results  presented  in  Reference  3. 
As  is  obvious  from  Table  10,  the  values  obtained  by  RRA  are  constant 
which  indicates  a  constant  spectral  shape.   However,  the  results  obtained 
from  this  work  range  from  .85  at  100  meters  to  .96  at  700  meters  which 
indicates  a  significant  change  in  the  spectral  shape  with  distance  as  is 
observed  in  Fig.  30. 
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Table  10.   Comparison  of  f  Values  for  the  8"  Shielded 
Source  Configuration 


Radial  Distance 
(meters) 

This 
Work 

RRA3 

100 

.845 

.89 

200 

.884 

.89 

300 

.914 

.89 

400 

.912 

.89 

500 

.922 

.89 

600 

.936 

.89 

700 

.958 

.89 

Finally,  Table  11  compares  the  f  values  obtained  for  the  16" 
shielded  source  configuration.   The  larger  variation  of  the  f  values 
for  the  shielded  source  configurations  is  easily  understood  from  Figs. 
38  and  39.   Since  the  shielded  source  configuration  spectra  fall  off 
almost  exponentially,  the  middle  and  high  energy  portions  of  the  spectra 
contribute  very  little  to  the  total  exposure  rate  thus  practically 
eliminating  the  mitigating  effect  on  the  upward  trend  of  the  f  values 
which  was  discussed  in  the  previous  section. 


Table  11.   Comparison  of  f  Values  for  the  16"  Shielded 
Source  Configuration 


Radial  Distance 
(meters) 

This 
Work 

3 
RRA 

100 

.855 

.885 

200 

.883 

.885 

300 

.908 

.885 

400 

.912 

.885 

500 

.918 

- 
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Again,  the  values  obtained  by  RRA  are  constant  with  source/ 
detector  distance  while  the  values  obtained  in  this  work  range 
from  0.86  to  0.92  over  the  500  meter  measurement  baseline.   Con- 
sidering the  slight  changes  in  the  spectral  shape  observed  in  Fig.  31, 
the  upward  trend  of  the  results  obtained  here  seem  well-justified. 
As  was  seen  in  Figs.  30  and  31,  the  number  of  photons  lying  immediately 
above  the  75  keV  maximum  decreases  with  increasing  source/detector 
distance.   Consequently,  the  amount  by  which  the  HPIC  overestimates  the 
total  exposure  rate  decreases,  thereby  causing  f  to  approach  unity. 

4.2  Comparison  of  Measured  and  Calculated  Exposure  Rate  Values 

The  exposure  rates  for  the  collimated  detector  geometry  were 
calculated  from  the  unfolded  Nal  spectra  using  Eq.  (4.1-4).   These 
values  were  then  reduced  by  2%  to  take  into  account  the  non-ideal 
collimator  effectiveness  described  in  Chapter  2.   The  value  of  this 
correction  for  the  non-ideal  behavior  of  the  detector  collimator  was 
estimated  by  comparing  the  area  under  the  actual  collimator  traverse 
curve  (Fig.  11)  to  that  expected  in  the  ideal  case3. 

Exposure  rates  were  calculated  from  the  DOT  3.5  fluxes  tabulated 
in  Reference  3,  again  using  Eq.  (4.1-4).   Exposure  rate  values  ob- 
tained from  these  calculated  spectra  are  also  given  in  Reference  3 
for  4tt  geometry.   In  general,  the  values  presented  in  this  work  are 
approximately  3%  lower  than  those  calculated  by  RRA3,  the  only 
exception  being  the  16"  shield  configuration  results  which  were  in 
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much  closer  agreement.   These  minor  discrepancies  are  probably  due  to 
differences  in  the  material  constants  (e.g.,  u/p,  etc.)  used. 

A  small  error  was  made  for  the  8"  shielded  source  configuration 

3 
in  the  DOT  3.5  transport  calculations  .   Instead  of  the  proper  concrete 

source  shield  thickness  of  21.03  cm  (8.28"),  20.03  cm  (7.89")  was  used. 

In  an  attempt  to  correct  for  this  oversight,  a  two-point  exponential 

fit  was  done  between  the  42.83  cm  (16.86")  and  the  20.03  cm  exposure 

rate  results  and  used  to  interpolate  a  value  for  a  21.03  cm  shield. 

For  those  source-detector  distances  where  the  42.83  cm  results  were 

not  available,  the  unshielded  configuration  results  were  used.  This 

correction  resulted  in  an  approximately  15%  reduction  of  the  20.03  cm 

values  when  the  interpolation  was  made  between  the  two  shielded  data 

points.   When  the  unshielded  and  20.03  cm  data  points  were  used,  the 

reduction  was  approximately  20%.   Those  values  obtained  using  the 

shielded  data  points  are  expected  to  be  more  accurate  since  the  source 

configurations  are  more  closely  related. 

4.2.1  Comparison  of  Collimated  Detector  Results 
A  comparison  of  the  exposure  rates  for  a  collimated  detector  are 
shown  in  Fig.  40  for  the  seven  source/detector  distances  used  in  this 
experiment.   The  solid  line  shows  the  trend  of  the  10  group  results 
while  the  dotted  line  follows  the  39  group  predictions.   The  quantity 
plotted  on  the  ordinate  is  the  exposure  rate  multiplied  by  the 
square  of  the  source/detector  distance,  normalized  to  a  unit  source 


Air  Mass  Thickness  (g/cm  ) 

Fig.  40.   Comparison  of  measured  and  calculated  exposure  rates 
for  the  collimated  detector  geometry. 
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2 
strength.   Multiplication  of  the  exposure  rate  by  r  removes  the  effect 

2 
of  the  geometric  attenuation  factor  (i.e.,  1/r  )  and  thus  the  plot  of 

the  results  is  a  fairly  slowly  varying  function. 

Since  the  air  density  was  not  constant  for  each  measurement, 

a  more  accurate  indication  of  the  distance  from  the  source  in  terms  of 

mass  attenuation  is  the  air  density,  p,  times  the  source/detector 

distance,  referred  to  as  the  air  mass  thickness.   The  density  of  the 

moist  air,  p  in  g/1,  with  relative  humidity,  R,  in  percent  and  at 

n  3 

temperature,  T,  in  K  was  calculated  as 

.  1.2929  (273.13)  r   P    0.3783es   fg  ] ,  , 

air         2         '1013  ll.i-ij 

where  P  is  the  atmospheric  pressure  in  millibars  and  e  is  the  saturation 

s 

27 
vapor  pressure  of  water  at  temperature  T  . 

As  can  be  seen,  the  10  group  results  underpredict  the  collimated 

geometry  exposure  rate  by  approximately  15  to  20%  for  the  unshielded 

and  8"  shielded  source  configurations.   The  errors  in  the  16"  shield 

configuration  results  are  closer  to  30%.   This  is  not  surprising  since 

as  was  discussed  in  Chapter  3,  the  calculated  10  group  energy  spectra 

fell  slightly  below  the  unfolded  Nal  results.   The  39  group  results 

agree  very  well  with  the  values  obtained  from  the  measured  Nal  spectra 

except  at  some  of  the  closer  measurement  locations  where  the  effect  of 

the  in-silo  scattered  component  is  significant. 
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One  source  of  error  in  the  10  group  DOT  results  Is  the  fact  that 
the  ground  was  not  considered  in  these  calculations.   Instead,  the  lower 
boundary  was  treated  as  a  vacuum.   It  has  been  estimated  that  the 

exclusion  of  the  ground  from  these  calculations  causes  an  under- 

3 
prediction  of  approximately  10%  . 

4.2.2  Comparison  of  4tt  Geometry  Results 

As  an  independent  verification  of  the  HPIC  measurements,  an  attempt 
was  made  to  infer  the  4u  geometry  exposure  rate  from  the  unfoided  Nal 
spectra.  This  technique  relies  heavily  upon  the  DOT  results  to  calcu- 
late a  geometry  correction  factor,  C,  which  is  defined  as 

c  ,  X  (collimated  geometry)  ^  2 

X  (4jt  geometry) 

Therefore,  the  exposure  rate  that  would  have  been  measured  by  the 

Nal  detector  in  4ir  geometry,  XT'1'.,  is  given  by 

;coll. 
•4it    -  Nal 
^al     C  (4-2-3) 

where  X^    is  the  exposure  rate  calculated  from  the  unfolded  Nal 
spectra  (i.e.,  collimated  geometry).   In  general,  C  is  a  function  of 
source/detector  distance  and  source  configuration. 

The  calculation  of  C  is  very  sensitive  to  any  errors  in  normali- 
zation of  the  DOT  calculated  exposure  rates.   If  the  error  in  the  magni- 
tude of  the  results  is  identical  for  the  collimated  detector  geometry 
results  and  the  4ir  values,  then  the  error  in  the  value  of  C  obtained  by 
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Caking  the  ratio  of  these  two  quantities  will  be  minimal.  However,  if  the 
normalization  error  is  different  for  the  two  calculated  exposure  rate 
values  then  the  error  incurred  in  the  calculated  value  of  C  will  be 
significant. 

Three  sets  of  C  values  were  calculated  in  this  analysis,  one  for 
each  source  configuration.   The  C's  used  in  the  inference  of  the 
unshielded  4tt  exposure  rates  were  calculated  at  100  meter  intervals 
by  taking  the  ratio  of  the  DOT  calculated  exposure  rates  for  the  colli- 
mated  and  4  it  geometry.  The  values  of  C  used  in  the  shielded  config- 
urations were  calculated  using  DOT  exposure  rates  for  the  shielded  results. 
The  results  of  these  calculations  are  shown  in  Table  12.  Quadratic  inter- 
polation was  used  between  the  values  tabulated  in  Table  12  to  obtain  values 
for  those  source/detector  distances  at  which  measurements  were  taken. 

Table  12.   Calculated  Geometry  Correction  Factors 
Shield  Thickne s s 


Radial  Distance 
(meters) 

0.0 

100 

.559 

200 

.523 

300 

.505 

400 

.496 

500 

.486 

600 

.485 

700 

.480 

16" 


.472  .499 

.455  .496 

.467  .454 

.461  .445 

.374  .424 
.347 
.313 


The  calculated  and  measured  4tt  geometry  exposure  rates  are 
shown  in  Fig.  41.  Again,  the  solid  line  follows  the  trend  of  the  10 
group  results  while  the  dotted  line  follows  the  39  group  predictions. 
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Comparison  of   measured   and   calculated   At  exposure   rates. 
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The  Reuter-Stokes  data  represents  an  average  of  at  least  two  and  at 
most  six  values.  As  in  the  collimated  geometry  results,  the  10  group 
values  underestimate  the  Reuter-Stokes  HPIC  values  by  approximately  10%  at 
all  locations  for  the  unshielded  configuration.  This  trend  also  occurs 
in  the  16"  shielded  configuration  results  where  the  differences  range 
from  20  to  25%  at  100  and  400  meters,  respectively. 

The  39  group  results  start  slightly  lower  than  the  10  group  values. 
An  error  of  almost  20%  is  observed  at  100  meters  when  the  39  group 
result  is  compared  to  the  Reuter-Stokes  datum.   At  500  meters  the  39 

group  value  is  approximately  4%  higher  than  the  measured  value. 

The  trend  observed  in  the  8"  shield  case  is  somewhat  different. 
Although  the  DOT  calculations  still  slightly  underestimate  the  inferred 
Nal  values,  the  Reuter-Stokes  results  are  much  lower  than  expected  at 
600  and  700  meters.  A  possible  explanation  for  the  sudden  decrease  in  the 
measured  exposure  rate  values  is  the  terrain  at  the  particular  measure- 
ment positions.   At  both  locations,  the  detector  was  positioned  on  a 
small  ridge.  At  the  600  meter  location,  small  valleys  were  situated 
on  either  side  of  the  detector  while  at  the  700  meter  location,  a  rather 
large  valley  was  positioned  directly  in  front  of  the  collimator  opening. 
The  absence  of  the  earth-air  interface  in  the  vicinity  of  the  detector 
resulted  in  a  significant  loss  of  ground  scattered  photons,  and  there- 
fore  an  approximately  10%  reduction  in  the  exposure  rate. 

The  agreement  between  the  Reuter-Stokes  HPIC  results  and  those 
values  inferred  from  the  unfolded  Nal  spectra  is  excellent.   In  most 
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cases  the  results  are  practically  Indistinguishable  from  each  other. 
The  only  instances  where  there  are  significant  differences  between  the 
two  sets  of  values  occur  at  600  and  700  meters  for  the  8"  shielded 
source  configuration.  Here  it  is  observed  that  the  Nal  results  con- 
tinue the  trend  established  by  the  values  at  the  closer  source/detector 
distances  lending  credibility  to  the  estimated  affect  of  ground- 
scattered  photons  on  the  total  exposure  rate.  The  angular  distribution 
of  the  low  energy,  ground  scattered  photons  could  be  such  that  the 
majority  of  them  are  shielded  from  the  Nal  detector  and  play  an  insignifi- 
cant role  in  the  collimated  geometry  results.   The  close  agreement 
between  these  two  sets  of  data  makes  a  clear  argument  for  the  accuracy 
of  the  experimental  results  and  the  deficiency  of  the  approximations 
used  in  the  photon  transport  calculations. 

The  almost  constant  separation  between  the  measured  and  calculated 
exposure  rate  values  causes  one  to  question  the  source  strength  values 
used  in  the  normalization  of  the  results  to  a  unit  source  strength. 
Figure  42  presents  the  4ir  exposure  rates  obtained  for  the  unshielded 
source  configuration.   Contained  in  Fig.  42  are  the  data  shown  in  the 
top  curve  of  Fig.  41,  but  now  the  sources  used  for  the  Nal  measure- 
ments are  identified.  A  comparison  between  the  results  obtained  with 

the  small  source  (10  Ci)  and  the  medium  source  (250  Ci)  is  observable 

2 
at  300  meters  (pr  =  34.32  g/cm  ).   Although  there  is  a  slight  discrepancy 

between  the  values,  it  is  certainly  no  larger  than  the  uncertainty  in 

the  reproducibility  of  the  results.   The  same  comments  are  applicable 

to  the  comparison  between  the  medium  and  big  (3800  Ci)  sources  at 
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2 
600  meters  (pr  =  68  g/cm  ).   Thus,  the  discrepancies  between  the  calculated 

and  measured  results  apparently  reflect  a  deficiency  in  the  calculational 

method  and  not  a  normalization  error. 

The  4tt  exposure  rate  values  for  the  unshielded,  8",  and  16"  shielded 

source  configurations  are  tabulated  in  Tables  13,  14,  and  IS,  respectively. 
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Table  13A.   Ait  Geometry  Nal  Exposure  Rates  for  the 

Unshielded  Source  Configuration 

2  • 

Air  Mass  Thickness                    r  -X 

(g/cm2)                     (m2-pR/h-Ci) 

11.375                       1.020  (5)* 

11.340                       9.859  (4) 

22.001                       1.047  (5) 

33.754                       7.676  (4) 

33.7A0                       7.536  (A) 

46.22A                         A. 868  (A) 

58.873                       2.786  (A) 

67.620                       1.889  (A) 

67.207                       2.032  (A) 

75.987                       1.377  (A) 

*1.020  (5)  =  1.020  x  105 
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Table  13B.   4ir  Geometry  DOT  Calculated  Exposure 
Rates  for  the  Unshielded  Source  Con- 
figuration 


Air  Mass  Thickness 

2 
r 

X 

(g/cm  ) 

On2- 

uR/h-Ci) 

(39  Group) 

(10  Group) 

11.44 

7.668 

(4) 

8.216  (4) 

22.88 

8.217 

(4) 

7.411  (4) 

34.32 

6.462 

(4) 

5.279  (4) 

45.76 

4.447 

(4) 

3.867  (4) 

57.20 

2.845 

(4) 

2.458  (4) 

68.64 

1.738 

(4) 

1.470  (4) 

80.08 

1.031 

(4) 

8.656  (3) 

7.668  (4)  =  7.668  x  10 


m 

Table  13C.  4tt  Geometry  Corrected 

HPIC  Exposure  Rate  Measure- 

ment»  i 

or  the 

Unshielded  Source 

Configuration 

Air  Mass  Thickne 

3S 

2 

r 

1 

(g/cm  ) 

2 
(m  -wR/h-Ci) 

11.58 

9.282 

A 

(4) 

11.44 

9.298 

(4) 

10.88 

9.123 

(4) 

10.88 

9.293 

(4) 

10.96 

9.663 

(4) 

11.20 

9.965 

(4) 

22.88 

9.788 

(4) 

21.92 

9.697 

32.88 

6.675 

(4) 

34.32 

9.957 

(4) 

34.32 

6.567 

(4) 

32.88 

6.833 

(4) 

45.76 

4.392 

(4) 

43.84 

4.947 

(4) 

57.20 

2.565 

(4) 

54.80 

2.917 

(4) 

67.44 

1.724 

(4) 

67.44 

1.824 

(4) 

65.76 

1.949 

(4) 

77.77 

1.200 

(4) 

76.72 

1.194 

(4) 

*9.282  (4)  -  9.282  x 

104 
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Table  14A.   4t  Geometry 

Nal 

Exposure  Rates  for  the  8" 

Shielded  Source 

Confi 

guration 

Air  Mass  Thickness 

r2-X 

(g/cm  ) 

(m2-yR/h-Ci) 

11.46 

7.923  (3)* 

11.57 

8.249  (3) 

11.68 
23.71 

7.696  (3) 
6.007  (3) 

35.03 

3.729  (3) 

34.58 

3.444  (3) 

35.28 

3.700  (3) 

45.73 

2.175  (3) 

57.67 

1.263  (3) 

68.68 

6.994  (2) 

68.60 

7.585  (2) 

78.59 

4.713  (2) 

*7.823  (3)  =  7.823  x  103 

Table  14B.  4ti  Geometry  DOT 

Calci 

jlated 

Exposure  Rates  for  the 

8"  Shielded  Sour 

ce  Configuration 

Air  Mass  Thickness 

r2-X 

(g/cm") 

2 
(m  •viR/h-Ci) 

11.44 

6.092  (3)* 

22.88 

4.553  (3) 

34 .  32 

2.826  (3) 

.'(5.76 

1.611  (3) 

52.70 

1.052  (3) 

68.64 

5.810  (2) 

80.08 

3.334  (2) 

*6.092  (3)  »  6.092  x  103 
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Table  14C.   Corrected  HPIC  Exposure  Rate  Measurements  for 

the 

8"  Shielded  Source  Configuration 

2 
Air  Mass  Thickness                 r 

X 

2                      2 
g(g/cm  )                  (m  -uR/h-Cl) 

11.08                      6.853 

(3)* 

11.62                     7.015 

(3) 

12.13                     7.350 

(3) 

12.13                     7.752 

(3) 

23.92                     5.778 

(3) 

22.24                      5.798 

(3) 

23.14                     5.793 

(3) 

35.74                     3.436 

(3) 

35.74                     3.536 

(3) 

33.49                     3.637 

(3) 

34.57                     3.686 

(3) 

47.64                     1.889 

(3) 

44.64                     2.321 

(3) 

46.08                     2.012 

(3) 

59.31                     9.588 

(2) 

55.63                     1.147 

(3) 

57.20                      1.053 

(3) 

66.51                     5.970 

(2) 

69.15                     5.379 

(2) 

78.14                     3.100 

(2) 

80.35                     3.120 

(2) 

*6.853  (3)  -  6.853  x  103 
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Table  15A.   4tt  Geometry 

Nal 

Exposure 

Rates  for  the  16"  Shielded 

Source  Configuration 

Air  Mass  Thickness 

r2-X 

(g/cm  ) 

(m  •uR/h-Ci) 

10.99 

7.136  (2)* 

11.95 

6.329  (2) 

11.96 

6.926  (2) 

11.06 

6.451  (2) 

23.67 

4.774  (2) 

35.09 

2.989  (2) 

34.83 

2.816  (2) 

46.34 

1.599  (2) 

56.74 

7.806  (1) 

56.71 

7.701  (1) 

57.21 

7.751  (1) 

*7.136  (2)  -  7.136  x  102 

Table  15B.  4ir  Geometry  DOT  Calculated 

Exp 

osure  Rates  for  the  16" 

Shielded  Source  Conf 

.gurat 

Lon 

Air  Mass  Thickness 

r2-X 

2 
(g/cm  ) 

2 
(m  -pR/h-d) 

11.44 

4.561  (2)* 

22.88 

3.162  (2) 

34.32 

2.007  (2) 

45.76 

1.077  (2) 

57.20 

5.638  (1) 

*4.561  (2)  -  4.561  x  102 
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Table  15C.   4n  Geometry  Corrected  HPIC  Exposure 

Rate  Measure- 

ments  for  the  16"  Shielded  Source  Configuration 

2 
Air  Mass  Thickness                   r 

X 

(g/cm2)                     (m2-uR/h-Ci) 

. 

11.10                     5.807 

(2)* 

11.18                       5.544 

(2) 

11.42                       5.588 

(2) 

22.19                     4.702 

(2) 

22.35                     4.621 

(2) 

23.01                     4.669 

(2) 

33.29                     2.732 

(2) 

33.64                     2.583 

(2) 

34.61                     2.523 

(2) 

44.52                     1.589 

(2) 

44.83                     1.473 

(2) 

46.15                     1.281 

(2) 

*5.807  (2)  =  5.807  x  102 

144 


5.0  Conclusions  and  Suggestions  for  Further  Study 

The  data  collected  in  this  experiment  represent  a  significant 
contribution  to  the  body  of  knowledge  in  the  area  of  radiation  shielding. 
These  data  should  prove  invaluable  in  the  development  of  new  methods  and 
the  verification  of  old  techniques  for  predicting  skyshine  exposure  rates. 
Furthermore,  the  results  of  this  experiment  may  be  used  directly  as 
design  data  in  the  construction  of  nuclear  facilities  where  skyshine 
radiation  is  of  concern. 

The  SEGO  unfolding  code  used  in  the  analysis  of  the  experimental 
data  appeared  to  work  very  well  on  the  smooth  spectra  encountered. 
The  comparatively  small  number  of  experimentally  determined  response 
functions  necessary  to  define  the  Hyodo  response  surface  simplified  the 
implementation  and  use  of  this  technique.   The  unfolded  results  also 
seemed  fairly  insensitive  to  minor  changes  in  the  response  matrix.   The  only 
problem  encountered  during  the  spectral  unfolding  was  the  iodine- 
escape  peak  which  had  to  be  artifically  suppressed  since  the  response 
surface  used  did  not  include  this  feature. 

The  agreement  between  the  unfolded  Nal  spectra  and  the  calculated 
spectra  was,  in  general,  very  good.  Excellent  results  were  obtained 
from  the  39  group  calculations  both  with  regard  to  spectral  shape  and 
normalization.   However,  the  10  group  calculations  were  somewhat  deficient 
in  describing  the  spectral  shape  and  as  a  consequence  underpredicted 
the  total  exposure  rate  by  approximately  20  to  30%. 
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The  excellent  agreement  between  the  corrected  Eeuter-Stokes 
high  pressure  ionization  chamber  (HPIC)  measurements  and  the 
collapsed  Nal  spectra,  both  independent  determinations  of  the  total 
exposure  rate,  lends  credibility  both  to  the  approximations  used  in 
the  calculation  of  the  energy  dependence  correction  factor,  f,  and 
to  the  accuracy  of  the  experimental  measurements  themselves. 

Several  portions  of  this  experiment  certainly  require  additional 
analysis.   Certainly  more  refined  transport  calculations  should  be 
performed  in  an  effort  to  resolve  the  apparent  discrepancies  between 
the  calculated  coarse  group  results  and  the  measured  values.   The 
effect  of  the  vacuum  boundary  condition  used  as  an  approximation  for  the 
actual  ground  interface  in  these  transport  calculations  should  be 
studied  in  some  detail.  Additional  verification  of  the  f  calculational 
procedures  used  in  this  study  could  be  obtained  by  performing  the  39 
group  calculations  for  the  shielded  source  configurations. 

Further  analysis  of  the  spectral  unfolding  technique  would  un- 
doubtedly prove  interesting.  Verification  of  the  unfolded  results 
could  possibly  be  obtained  through  the  use  of  a  different  unfolding 
procedure.   There  also  appears  to  be  a  slight  philosophical  weakness 
in  the  underlying  assumptions  of  the  SEGO  code.   As  formulated  in 
Eq.  (3.1-6),  the  system  resolution  is  a  function  of  the  incident  photon 
energy  alone  while  physically  the  smoothing  is  a  function  of  the  amount 
of  energy  actually  deposited  in  the  crystal.  Thus  Eq.  (3.1-6)  should 
read 
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E 

t  max 

8^(1')  =    dE   dE"  U(E',E")  S(E",E)  i=l,2,...N  .  (5.0-1) 

AEJ    0 

l 

An  effort  should  be  made  to  reconcile  the  two  approaches. 

Finally,  much  work  remains  in  the  development  of  accurate  point 
kernel  techniques  to  predict  sky shine  exposure  rates.  The  use  of  a 
buildup  factor  in  a  finite  media  problem  will  inevitably  lead  to  some 
error  as  the  source/detector  distance  increases.  Perhaps  a  more 
appropriate  buildup  factor  could  be  constructed  from  the  available 
transport  theory  results. 
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Abstract 

A  benchmark  gamma-ray  skyshine  experiment  was  performed  to  obtain 
skyshine  exposure  rate  and  spectral  data  against  which  numerical  models 
may  be  compared.   To  obtain  skyshine  radiation  fields  of  measureable 
intensity  over  the  700  meter  measurement  baseline,  three   Co  sources 
were  used  of  nominal  strengths  10,  250,  and  5000  Ci.   Three  source 
configurations  were  employed;  first  with  the  source  radiation  collimated 
into  a  vertical  conical  beam  and  then  with  the  sources  covered  with 
eight  and  sixteen  inches  of  concrete. 

Spectral  measurements  were  made  at  100  meter  intervals  from  the 
Co-60  source  with  a  collimated  Nal(Tl)  scintillation  spectrometer.  The 
detector  was  housed  in  an  eight  ton  collimator-shield  assembly  and 
mounted  on  a  semi-trailer.   Total  skyshine  expsosure  rates  were  also 
measured  over  the  700  meter  baseline  with  a  ten- inch  diameter,  25  atmo- 
sphere argon  filled  Reuter-Stokes  high  pressure  ionization  chamber. 

The  measured  Nal  spectra  were  unfolded  with  the  SEGO  code  using 
experimentally  determined  response  functions.   These  unfolded  spectral 
results  are  compared  with  results  of  detailed  DOT-3.5  transport  calculations 
and  found  to  agree  to  within  20%. 

A  factor  was  calculated  using  the  unfolded  Nal  spectra  to  correct 
the  Reuter-Stokes  detector  measurements  for  its  non-ideal  energy  response. 
Additionally,  the  unfolded  Nal  spectra  are  used  to  obtain  an  independent 
determination  of  the  total  skyshine  exposure  rate.   These  independently 


measured  values  agreed  to  within  5-10%,  while  values  calculated  by 
the  DOT  code  were  consistently  10  to  20%  lower. 


